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ABSTRACT 

Cepheus B (Cep B) molecular cloud and a portion of the nearby Cep 0B3b 
OB association, one of the most active regions of star formation within 1 kpc, has 
been observed with the ACIS detector on board the Chandra X-ray Observatory. 
We detect 431 X-ray sources, of which 89% are confidently identified as clus- 
tered pre-main sequence (PMS) stars. Two main results are obtained. First, we 
provide the best census to date for the stellar population of the region. We iden- 
tify many members of two rich stellar clusters: the lightly obscured Cep 0B3b 
association, and the deeply embedded cluster in Cep B whose existence was pre- 
viously traced only by a handful of radio sources and T Tauri stars. Second, we 
find a discrepancy between the X-ray Luminosity Functions (XLFs) of the Cep 
0B3b and the Orion Nebula Cluster (ONC). This may be due to different Initial 
Mass Functions (IMFs) of two regions (excess of ~ 0.3 Mq stars), or different 
age distributions. 

Several other results are obtained. A diffuse X-ray component seen in the 
field is attributed to the integrated emission of unresolved low mass PMS stars. 
The X-ray emission from HD 217086 (07n), the principle ionizing source of the 
region, follows the standard model involving many small shocks in an unmag- 
netized radiatively accelerated wind. The X-ray source #294 joins a number of 
similar superflare PMS stars where long magnetic structures may connect the 
protoplanetary disk to the stellar surface. 
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1. Introduction 

Ccphcus B (Ccp B), the hottest CO component of the Cepheus molecular cloud (Sargent 
1977, 1979), is located near the northwestern edge of the cloud at a distance of 725 pc from 
the Sun (Blaauw et al. 1959; Crawford & Barnes 1970). Outside Cep B, to the north and 
west of the cloud, lies the stellar Cep 0B3 association, composed of two subgroups of stars 
of different ages with the youngest (Cep 0B3b) lying closer to the molecular cloud (Blaauw 
1964; Jordi et al. 1996) (see Figure 1). 

Sargent (1979) considered Cep 0B3 to be a good example of sequential star formation, 
following the model of Elmegreen & Lada (1977). The interface between the molecular cloud 
and the OB star association is clearly dehneated by the optically visible HI I region Sharpless 
155 (S 155), where cloud material is ionized and heated by the radiation field of the youngest 
generation of the Cep OB association (Panagia & Thum 1981). The photodissociation region 
(PDR) at S 155 is favorably oriented to reveal the progression of star formation. As the 
surface of the cloud is being eroded by the early-type stars, the cloud edge moves eastward 
across the observer's field of view with new stars emerging from the obscuring molecular 
cloud. Radio continuum (Felli et al. 1978; Testi et al. 1995), far-infrared, CO (Minchin et 
al. 1992), and near- infrared studies (Moreno-Corral et al. 1993; Testi et al. 1995) reveal a 
few young stars embedded in the Cep B cloud behind the PDR. Sources located on the edge 
of the molecular cloud and exhibiting high extinction were suggested to represent a third 
generation of star formation triggered by the HII region shock propagating into the cloud. 

The low mass pre-main sequence (PMS) stellar population of Cep 0B3b can be ef- 
fectively discriminated using X-ray emission due to the heavy contamination in optical and 
near-infrared (ONIR) surveys by both foreground and background Galactic field stars. X-ray 
surveys are complementary to ONIR surveys because they trace magnetic activity (mainly 
plasma heated in violent magnetic reconnection fiares) rather than photospheric or circum- 
stellar disk blackbody emission, and PMS X-ray emission is elevated 10^ — 10'^ above main 
sequence (MS) levels (Preibisch et al. 2005). A ROSAT X-iaj study of Cep 0B3b identi- 
fied over 50 likely PMS stars with Ay < 3 mag and Lx ^ 10^° ergs s^^ (Naylor & Fabian 
1999). Ten of these ROSAT sources are coincident with optical sources, and four of these 
are spectroscopically confirmed classical and weak-lined T Tauri stars (Pozzo et al. 2003). 

Much more can be achieved with the current generation of X-ray telescopes, particularly 
the Chandra X-ray Observatory with its excellent high-resolution mirrors and low-noise de- 
tectors. Studies of the Orion Nebula region demonstrate that Chandra images can penetrate 
up to Ay ^ 500 mag into the cloud (Getman et al. 2005; Grosso et al. 2005), considerably 
deeper than near-infrared surveys in the JHKL bands. Chandra is also effective in resolving 
crowded fields down to ~ 0.7" scales (Getman et al. 2005b). In X-rays, OB stars are often 
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not much brighter than PMS stars, so multiple systems associated with OB stars can be 
identified (Stelzer et al. 2005). Chandra X-ray studies are particularly efi^ective in uncov- 
ering heavily obscured low-mass cloud populations and in discriminating PMS populations 
Irom unrelated older stars. 

We present here a Chandra study of the region using its Advanced CCD Imaging Spec- 
trometer (ACIS) detector. Over 380 PMS stars in two rich clusters arc found, concentrated 
in the lightly obscured Cep 0B3b cluster and a previously unknown embedded chister, and 
we infer the existence of an additional 700 — 800 low-mass members. Only ~ 7 — 9% of the 
detected X-ray sources arc attributable to contaminating Galactic and cxtragalactic popu- 
lations. The observation and source lists are described in §2. Identification with 2MASS 
counterparts and their NIR properties are considered in §3. X-ray fluxes and their errors are 
determined in §4. X-ray sources unrelated to the cloud are discussed in §5. The discovery 
of two rich clusters is described in §6. Their stellar populations are discussed in §7 using 
measured X-ray luminosity functions (XLFs) and Initial Mass Functions (IMFs). The origin 
of the soft diffuse X-ray emission discovered here is discussed in §8. Finally, we address 
two individual stars of particular interest: HD 217086, the most massive star in Cep 0B3b 
(§9.1), and Chandra source #294 which exhibited an unusually luminous and hot X-ray flare 
(§9-2). 



2. The Chandra observation and source list 

2.1. Observation and data reduction 

The observation of Cep B and Cep 0B3b was obtained on 2003 March 11.51 — 11.88 
with the ACIS camera (Garmire et al. 2003) on-board Chandra (Weisskopf et al. 2002). We 
consider here only results arising from the imaging array (ACIS-I) of four abutted 1024 x 1024 
pixel front-side illuminated charge-coupled devices (CCDs) covering about 17' x 17' on the 
sky. The S2 and S3 detectors in the spectroscopic array (ACIS-S) were also operational, but 
as the telescope point spread function (PSF) is considerably degraded far off-axis, the S2-S3 
data are omitted from the analysis. The aim point of the array is 22'^56™49'!4, +62°39'55'.'6 
(J2000), and the satellite roll angle (i.e., orientation of the CCD array relative to the north- 
south direction) is 7?9. The total net exposure time of our observation is 30 ks with no 
background flaring or data losses. 

Data reduction, source detection, and extraction have been done using procedures sim- 
ilar to those described in detail by Gctman et al. (2005). Briefly, the data were partially 
corrected for CCD charge transfer inefficiency (Townsley et al. 2002), cleaned with "grade". 
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"status", and "good-time interval" filters, trimmed of background events outside of the 
0.5 — 8.0 keV band, and cleaned of bad pixel columns with energies of < 0.7 keV, left by the 
standard processing (see Appendix B of Townsley et al. (2003) for more details). Based on 
several dozen matches between bright Chandra and 2MASS point sources a correction to 
the absolute astrometry was applied and the slight PSF broadening from the Chandra X-ray 
Center's (CXC's) software randomization of positions was removed. Source searching was 
performed with data image and exposure maps constructed at three spatial resolutions (0.5, 
1.0, and 1.4" per pixel) using the CI AO wavdetect tool. Those procedures were performed 
with CIAO software package 3.1, LHEASOFT 5.3, the Penn State CTI corrector version 
1.42, and the ACISExtract (AE) package version 3.67. The latter two tools were developed 
at Penn State^. 



2.2. Source detection and characterization 

Figure 2 shows the resulting smoothed image of the ACIS-I field. More than 400 point 
sources are seen^. A wavelet-based source detection procedure (wavdetect) is applied. It is 
followed by visual examination to locate other candidate sources, mainly close doubles and 
sources near the detection threshold. Using AE, photons are extracted from the polygonal 
contours of ~ 90% encircled energy using detailed position-dependent models of the PSF. 
Background is measured locally in source-free regions. The list of candidate sources is 
trimmed to omit sources with fewer than three background-subtracted counts in the PSF 
(i.e. NetFull/ PSF < 3; these and other quantities are defined in the notes of Table 1). 
This resulted in a catalog of 431 sources, which are hsted in Table 1. Twenty-seven of the 31 
weakest sources with 3 < NetFull/ PSF < 5 counts (median source significance Signif — 
1.1, which corresponds (Getman et al. 2005) to 2.2 a detection threshold of PWDetect) have 
near-infrared (NIR) counterparts (§3); we thus believe most of these arc real despite their 
marginal statistical significance. One hundred and eight sources with 5 < NetFull /PSF < 
10 have median Signif = 1.8 (3.6 cr detection threshold of PWDetect). 

The AE package provides a variety of source characteristics including celestial position, 
off-axis angle, net and background counts within the PSF-based extraction area, source sig- 



^ Descriptions and codes for CTI correction and ACIS Extract can be found at 
http://www.astro.psu.edu/users/townsley/cti/ and http://www.astro.psu.edu/xray/docs/TARA/ae_users_guide.html, 
respectively. 

^We do not show here the raw image of the field as it gives a misleading impression of high background 
and low resolution. 
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nificance assuming Poisson statistics, effective exposure (corrected for telescope vignetting, 
satellite dithering, and bad CCD pixels), median energy after background subtraction, a 
variability indicator extracted from the one-sided Kolmogorov-Smirnov statistics, and occa- 
sional anomalies related to chip gap or field edge positions. These are all reported in Table 
1; see Getman et al. (2005) for a fuller description of these quantities. 

3. 2MASS counterparts and their properties 

3.1. Identification with catalogued 2MASS sources 

An automated cross-correlation between the Chandra source positions and 2MASS NIR 
catalogs was made using a search radius of 1" within ~ 3.5' of the ACIS field center, and a 
search radius of 2" in the outer regions of the field where the Chandra PSF deteriorates. To 
evaluate the merits of this automated procedure, we performed a careful visual examination 
of each source. In 99% of the cases, the identifications are clear and unambiguous. Seventeen 
matches near the edges of the field with Chandra-2MASS separations greater than 2" were 
added to the catalog. The result is that 385 out of 431 (89%) Chandra sources have NIR 
counterparts in the 2MASS catalog (Table 2). 

After a correction to the Chandra boresight of AR.A. = -1-0.13" and ADecl. = —0.14", 

mentioned in section 2.1, to bring it into the 2MASS/Hipparcos reference frame, the median 
Chandra-2MASS separations arc better than 0.21" at the central part of the field (ofi'-axis 
angle < 3.5'), and 0.47" at larger ofi^-axis angles. This excellent astromctric agreement 
between 2MASS and Chandra positions is similar to that achieved in the Orion Nebula 
Cluster field (Getman et al. 2005). 

3.2. X-ray sources with un-catalogued 2MASS counterpcirts 

Of the 46 X-ray sources not associated with 2MASS catalogued stars 14 are soft {MedE < 
2 keV) and 32 are hard {MedE > 2 keV). Visual examination of the 2MASS atlases shows 
that 12 soft (## 78, 103, 170, 222, 242, 243, 249, 296, 309, 340, 343, 374) and 2 hard (## 
395 and 403) X-ray sources without identification lie very close to bright NIR stars. As an 
example, expanded views of the 2MASS atlas around 12 soft X-ray sources are shown in 
Figure 3. The nearby catalogued 2MASS stars are often fiagged as an extended NIR source. 
In the case of 242 and 243, Chandra found two un-catalogued companions to a bright 
2MASS star. X-ray source # 340 is suspiciously close to an unregistered bright 2MASS 
source, but does not coincide with it; this may be emission from an HH shock, but the lack 
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of optical data does not allow us to speculate on it. The only two soft X-ray sources that are 
truly not seen in NIR 2MASS arc 26 and 214. These two are the hardest among the soft 
sample and perhaps are too absorbed to be detected in NIR. Later in the membership study 
(§5) we consider these 14 soft sources, including 7^7^ 26, 214, 340, and two hard sources, #7,^ 
395 and 403, as likely members of the Cepheus populations. 

3.3. NIR properties of X-ray detected sources 

Figures 4 and 5 compare NIR properties of all cataloged 2MASS sources in the ACIS-I 
field with the 385 detected in X-rays (not including un-cataloged 2MASS sources from §3.2). 
The latter sources are clustered in NIR color space at Ay ~ 2 — 4 mag with only ~ 10 
having NIR colors consistent with foreground main sequence stars. NIR colors for ~ 30, 8% 
of the X-ray selected population, show i^-band excesses that indicate the presence of dusty 
protoplanetary disks. 

Pozzo et al. (2003) found the age of ~ 1 Myr for four out of six PMS candidates within 
their la and 2b fields (see their Figure 15), which partially overlap our ACIS-I field. In Figure 
5 we show the 1 Myr isochrone in the J vs. {J — H) color-magnitude diagram with reddening 
loci for PMS stars of different masses (see footnote 4 for references) . We compare locations 
of our X-ray sources on the color-magnitude diagram with the 1 Myr PMS isochrone and 
estimate their masses, bolometric luminosities, and local absorptions. Most of the X-ray 
sources are consistent with PMS stars with masses of 0.1 — 1 M© at a distance of 725 pc 
(Testi et al. 1995) subject to 1 < Ay < 3 mag absorption. Several dozen appear to be more 
obscured PMS stars with 5 < Ay < 15. 

Table 2 provides information on the 385 2MASS counterparts, their offsets from Chan- 
dra positions, photometric magnitudes, photometry quality and confusion flags. With the 
assumption that most of the detected X-ray sources are ~ 1 Myr old PMS stars we also give 
our estimates of stellar masses, absorptions, and bolometric luminosities based on the NIR 
color-magnitude diagram in Figure 5. 

3.4. Understanding the background stellair population 

We examined Monte Carlo simulations of the Galactic stellar population expected in 
the direction of our ACIS field {I, h) — (110.2, -1-2.7) based on the stellar population synthesis 
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model of Robin et al. (2003)^. Within the ACIS field-of-view, the Besangon model predicts 
~ 100 - 150 MS foreground {d < 725 pc) and ~ 1400 - 3100 background MS, giants, and 
sub- giants within the 2MASS sensitivity hmit around J < 16 — 17 mag. 

The simulations help us to understand the bimodal distribution of 2MASS stars in the 
color-magnitude diagram (this was previously noticed from V versus V — I optical diagrams 
by Pozzo et al. (2003)). The model indicates that, due to its sensitivity limits, 2MASS 
will preferentially detect the M, K, and G MS stars at distances of ~ 0.8 — 3.6 kpc, but will 
detect large populations of A and F MS stars, and sub-giants and giants of all spectral types, 
located at distances of ~ 5.4 — 5.7 kpc. Loci for these sub-populations arc shown in the NIR 
color-magnitude diagram of Figure 5^. The region occupied by most of our X-ray sources, 
PMS stars with Ay ~ 2 — 4 mag, is also the characteristic location for ~ 400 — 500 older 
(> 1 Gyr) background giants. None of those giants is expected to emit X-rays detectable in 
our observation (§5.3). 



^This was computed by the Web service provided by the Robin et al. (2003) group at Besangon at 
http://bison.obs-bcsancon.fr/modele. In addition to the smooth absorption component provided in the 
model, we added a cloud at d = 725 pc with Ay ~ 2.5 mag, as inferred from the average NIR colors of 
Chandra sources shown in Figures 4-5. As we have not considered the patchy denser molecular cloud which 
covers ~ 1/4 of ACIS-I FOV, the inferred background population may be somewhat overestimated. We refer 
to this hereafter as the Besangon model. 

^We give here details concerning the J vs. (J — H) color magnitude diagram in Figure 5. Black dots 
show Chandra sources with 2MASS counterparts, and grey dots show all 2MASS sources in the field. The 
ZAMS (thick black dotted line on left) is obtained from Tables 7.6 and 15.7 of Cox (2000). The IMyr 
isochrone (thick black line on right) is from Baraffc et al. (1998) for 0.02 M© ^ M ^ 1.4 M© and Siess 
et al. (2000) for IAMq ^ M ^ 7.0 M©. Reddening vectors of Ay ^ 10 mag (dashed lines) arc given for 
0.075 Mq, 0.1 M0, 0.2 Mq, 0.5 Mq, 1.0 Mq, 2.0 Mq, and 4.0 Mq. Overlayed on this color-magnitude diagram 
are outlines of the distributions of expected foreground and background populations obtained from Besangon 
simulations of the Galactic stellar populations in the direction of Cepheus. They clearly explain the two- 
horn shape of the color distribution of 2MASS stars. The horn around J — ~ 0.5 — 0.8 consists of G-K-M 
main sequence background stars (red locus), A-F main sequence background stars (green locus), background 
sub-giants (luminosity class IV, purple locus), younger (< 1 Gyr) background giants (luminosity class III, 
black dashed locus around J — H = 0.5), and foreground main sequence stars (blue locus). The horn aroimd 
J —H ~ 0.8— 1.2, which includes most of the Chandra sources, has older background giants (> 1 Gyr, black 
dashed locus around J — H = l.Q). The solid black locus covering heavily obscured stars with Ay > 5 mag 
contains a mixture of embedded PMS stars and highly reddened background main sequence stars, sub-giants 
and giants. 



-8- 



4. X-ray spectroscopy and fluxes 

The XLF of stars in young stellar clusters exhibits a range of 10^ in X-ray luminosity 
(Feigelson et al. 2005). In relatively short exposures of distant clusters, the majority of 
sources have very few counts while a few may have thousands of counts. Consequently, no 
single method is well-adapted to characterization of source spectra and fluxes. We describe 
here three different approaches applied to three flux levels, which are summarized in Figure 
6: for 174 sources with > 20 counts (sixth row of the Figure), we perform regressions to 
spectral models giving direct estimations of plasma temperatures, line-of-sight absorption, 
and broad-band fluxes; for 212 weak sources with absorptions estimated from NIR colors^ 
(seventh row of the Figure), we perform regressions with plasma temperature as a parameter; 
and for 45 weak X-ray sources without Ayj^ estimates (eighth row of the Figure), we obtain 
non-unique spectral fits for the sole purpose of estimating observed broad-band fluxes. This 
last category of sources mostly consists of sources unrelated to the Cepheus star formation 
region (§§5.1 and 5.2). Resulting spectral fits and broad-band luminosities appear in Table 
3. Following Getman et al. (2005), we report observed hard- band (2 — 8 keV, logL/j) and 
full-band (0.5 — 8 keV, logL^) luminosities, and intrinsic luminosities after correction for 
absorption (log L^^ci log -^t,c) ■ We omit soft band luminosities due to the absorption seen in 
virtually all sources. 

4.1. Spectral fits for sources with > 20 counts 

For 174 sources with greater than 20 net counts (154 hghtly absorbed and 20 heavily 
absorbed), automated spectral fitting was performed using AE and XSPEC (Arnaud 1996). 
We applied a one-temperature optically thin thermal plasma emission model (Smith et al. 
2001, APEC) assuming a uniform density and 0.3 times solar elemental abundances which 
is typical of young stellar X-ray sources (e.g. Imanishi, Koyama, & Tsuboi 2001; Feigelson 
et al. 2002a). To explore the whole range of various alternative models that successfully fit 
the data we performed spectral fits on grouped (with different minimum number of grouped 
counts, starting from 1, using the x^-statistic^) as well as on un-grouped (using the C- 
statistic) spectra. We then performed a visual inspection of all spectral fits to choose the 



^Nh,ir was derived from Avjr reported in Table 2 using the relationship Nh = 2.0 x 10^^ Ay (Ryter 
1996). 

^We tested x^-statistics with different error weighting methods, such as standard, gehrels, and churazov, 
and found that the best-fit values of the model parameters remain similar for a spectrum of the same 
grouping. 
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"most appropriate" spectral fit based on subjective judgment as outlined below. 

For lightly absorbed sources (roughly, sources with median energy MedE < 2 keV), 
data can often be successfully fit with two or more qualitatively different spectral models. 
Spectral fits of five sources with 20 — 23 net counts are shown in Figure 7 to exemplify this 
common ambiguity in spectral fitting. Sources can be often successfully fit by two different 
types of models: high and low /cT, or lower with higher kT . Source ^ 167 can be fit 
by three different classes of models. This problem was also noted by Getman et al. (2005) 
for the Orion Nebula Cluster (COUP project) where many more photons are available. In 
these situations, wc adopted a model where the plasma energy most closely matches the 
typical energy of PMS stars in Orion: for unabsorbed stars with NIR counterparts, the 
median plasma temperature (from one-temperature fits) is kT ~ 1.6 keV; for absorbed stars 
with NIR counterparts, kT ~ 2.6 keV; and for absorbed stars without NIR counterparts, 
kT ~ 4.4 keV (Getman et al. 2005b). In Figure 7, we thus preferred the right-hand panel 
fits. For highly absorbed sources (all except 4 with > 30 net counts), data with different 
groupings are successfully fit by similar class of models (Figure 8), and we arbitrarily report 
a fit from one of the available groupings. 



4.2. Spectral fits for sources with < 20 counts 

Figure 9 compares absorption column densities derived from our X-ray spectra of bright 
{NetFuU > 20 cts) sources with those derived from the 2MASS J vs.{J—H) color magnitude 
diagram. NIR sources were dereddened to the 1 Myr PMS isochrone. We see that the X- 
ray and NIR absorption columns are statistically equivalent; no bias is seen and scatter 
is generally within ±0.3 in log Nh''^. On this basis, we can reliably substitute log Nhjr 
for log NH,xray in spectral fitting of 212 fainter sources (177 soft and 35 hard) which have 
2MASS NIR photometry. 

To estimate plasma temperatures in these faint sources, we use the median energy of 
the extracted photons. The COUP study showed that MedE is a predictor of log Nh in 
strong PMS stars superposed by a weaker dependence on kT (Feigelson et al. 2005). To test 
this method, we simulated a grid of absorbed plasma models (wabs*APEC) using the fakeit 
command in XSPEC. We then passed simulated spectra through AE to perform photometric 
analysis including calculation of MedE. For the 212 faint sources with NIR photometry, we 
fixed column densities to Nh,ir and placed them on the simulated {Nh versus MedE) plane 
to roughly estimate their plasma energies. With the Nh and kT parameters frozen we then 



''A similar correlation exists for COUP data, but it has not been reported. 
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performed XSPEC fits to derive the spectrum normalization and to calculate broad-band 
X-ray fluxes. Examples of these spectral fits are shown in Figures 10 and 11. The resulting 
spectral fits appear satisfactory, even for very weak sources of as few as 3-5 net counts. 

Finally, we consider the 45 weak sources (17 soft and 28 hard) that lack sufficient 2MASS 
photometry to estimate Nh,ir- Most of those sources are probably unrelated to the Cepheus 

region (sec §5). Here we report one of the possible spectral fits from a wide range of models in 
order to integrate under the model and to estimate observed broad-band fluxes. The spectral 
models for these 45 weak sources are used here simply as spline flts with no astrophysical 
meaning, mainly to provide future researchers with the ability to derive much more reliable 
observed X-ray fluxes. Table 3 reports thermal models for all 45 sources; Table 4 reports 
power law models for 21 of these sources (see last row in Figure 6) that are AGN candidates 
(see §5.1). 

4.3. Flux uncertainties 

Formal errors on X-ray luminosities for young stellar members of the Cepheus region 
have been estimated through Monte Carlo simulations. For each source: (1) we randomly 
draw a distribution of log Nh, comprised of 1000 points, with the mean equal to that reported 
from spectral analysis and a ~ 0.3, obtained in §4.2; (2) we randomly draw a count rate 
{CR) Poisson distribution of 1000 points with the mean equal to the measured value; (3) 
using the MARX ray-trace simulator, we simulate a very bright source with spectral source's 
characteristics, N^, kT, we slice the simulated source on thousand of chunks, comprising 
the number of counts of the real source, and derive the characteristic distribution of source's 
median energy (MedE); (4) we pass the obtained distributions of logNu, CR, and MedE 
through the simulated database of optically thin plasma model (see §4.2) to report the 
standard deviations of X-ray fluxes in Table 3. 

The simulations indicate that 1 a errors on MedE are in general better than 0.2 kcV for 
brighter {NetFull > 20 counts) sources, but become as large as ~ 1 keV for weaker sources. 
X-ray observed luminosities Lt {L^, Lh^c) are known better than 0.2 (0.3,0.3) dex for brighter 
sources, but become as much uncertain as 0.3 (0.6,0.6) dex for weaker sources (Figure 12). 
The uncertainty on intrinsic total band luminosity Lf c in addition to count rate is affected 
by absorption, but in general for most of the sources it is better than 0.4 dex. 

The scientific reliability of these luminosities, however, is usually lower than the formal 
statistical uncertainty for several reasons. First, all PMS stars exhibit strong temporal 
variabihty. From Orion studies, we estimate that ~10% of the stars will be caught in a fiare 
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state during the 8 hours of Cepheus B exposure (Wolk et al. 2005). Second, the Lt^c will 
systematically underestimate the true emission for obscured or embedded populations, due 
to the absence of soft band spectral components. Third, if the distance of 725 pc adopted in 
this paper were revised to 800 pc, quoted in Pozzo ct al. (2003), then the PMS stellar X-ray 
luminosities should systematically be increased by ~ 0.1 dex. 

5. X-ray sources unrelated to the cloud 

5.1. Unidentified X-ray sources and extragalactic contamination 

Out of 83 "hard" {MedE > 2 keV, Figure 6) sources 30 have not been identified with 
2MASS NIR counterparts (this does not include ## 395, 403 from §3.2). Twenty four of 
these sources without 2MASS counterparts lie outside of the Cep B molecular cloud, which 
covers the southeastern ACTS chip (CCDO) (Figure 13(a)). We argue here these are most 
likely extragalactic sources. No more than 3 extragalactic sources are likely present behind 
the dark cloud itself; most of the hard sources in the cloud arc concentrated around the 
PDF ionization front (see Figure 13(a)) or in the new heavily absorbed cluster found in the 
ACIS image (§6). We therefore omit consideration of extragalactic contamination on the 
southeastern ACIS chip. 

We evaluate the expected contamination by extragalactic X-ray sources using the meth- 
ods of Getman et al. (2005b). Confining analysis to the hard (2 — 8) keV band, we construct 
Monte Carlo simulations of the extragalactic population by placing artificial sources ran- 
domly across the field. Incident fiuxcs are drawn from the X-ray background log N — log S 
distribution described by Moretti et al. (2003) assuming a power law photon index distribu- 
tion consistent with flux-dependencies described by Brandt et al. (2001). The spectrum of 
each simulated source was passed through a uniform absorbing column density log Nh ~ 22.2 
cm~^. This is the estimated sum of two components: absorption of log Nh ~ 21.7 cm~^ from 
Earth to the Cepheus region taken from the median of soft X-ray sources away from the dark 
cloud and a HI column density of log Nh ~ 22.1 cm~^ through the entire Galactic disk in 
the direction of Cepheus (Dickey & Lockman 1990)®. After applying local background levels 
found in our ACIS image, we then remove very weak extragalactic sources that would have 
fallen below our source detection threshold. 

The grey lines in Figure 14 show the predicted — log S distribution of extragalactic 



^This was obtained from NASA's HEASARC tool located at littp://heasarc.gsfc. nasa.gov/cgi- 
biii/Tools/w3nh/w3nh.pl 
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sources from several typical simulations with different detection thresholds. Here N is the 
total number of extragalactic sources expected to be detected and S = Fh is the observed 
source flux in the hard band. This prediction is compared with the observed hard fluxes of 
24 Cepheus hard unidentified sources (black solid line in Figure 14; values given in Tabic 
4). The simulations predict that a range of 7 — 17 extragalactic sources would have been 
detected in these three CCD chips by our data processing, which is roughly half of the 24 
observed unidentified hard Cepheus sources. The others may be physically associated with 
the Cepheus region — for example new embedded sources in the molecular cloud, which may 
extend towards the southwestern (part of CCDl) and northeastern (part of CCD2) parts of 
the ACIS-I field (Figure 13(a)). Since the observed source counts he somewhat below the 
simulations for 5* < 1 x 10~^^ erg s~^ cm~^, it is plausible that 3 — 5 of the hard sources with 
2MASS counterparts are bright extragalactic sources rather than PMS stars in the cloud. 

As we know that roughly half of the 24 hard unidentified sources are extragalactic 
contaminants, but we can not make individual determinations, we flag all 24 sources as 
possible AGN contaminants (flag 'a' in Table 3). They are omitted from further consideration 
of the Cepheus PMS population. 

5.2. Foreground stellar contamination 

Figure 15 shows an expanded view of a portion of the NIR color-magnitude diagram 
in Figure 5. Thirteen lightly- absorbed Chandra stars lie between the 1 Myr isochrone and 
the ZAMS within the locus of Besangon simulated foreground stars and thus are unlikely 
to be Cep B cloud members. These are also indicated schematically at the bottom left of 
the X-ray source classiflcation diagram shown in Figure 6. Although X-ray surveys of star 
forming regions are very effective in preferentially detecting PMS stars over MS stars, even 
a very low detection rate of older stars will lead to interlopers in the X-ray image due to 
the huge number of field stars. Field stars detected in X-ray surveys will preferentially be 
younger members of the disk population, as stellar X-ray emission decays rapidly after ~ 1 
Gyr (Feigelson et al. 2004; Preibisch & Feigelson 2005b). 

Half of the field stars in the Besangon simulation have ages greater than 2 Gyr while 
half have ages ranging from 0.15 to 2 Gyr. Most reside at distances from 400 to 725 pc; 
65% are M stars, 24% are K stars, 9% are G stars, and 2% are more massive. We convolved 
the Besongon model populations with the X-ray luminosity functions of stars in the solar 
neighborhood measured from ROSAT surveys (Schmitt et al. 1995; Schmitt 1997; Hiinsch 
et al. 1999). Luminosities were adjusted to account for the different ROSAT a,m\ Chandra 
spectral bands following the stellar hardness-luminosity relation (Giidel et al. 1998). The 
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resulting predicted X-ray flux of tlie foreground Besangon sample in the (0.5 — 8.0) keV band 

is Fx ~ 2 ± 1 X 10~^^ erg cm~^. Following the same procedure we used for extragalactic 
sources (§5.1), wc applied our ACIS detection process to these simulated field stars. The 
results suggest that 5—10 field stars may be detected in the AClS-1 field. A typical Monte 
Carlo run predicts 2 F stars, 3 G-K stars, and 4 M stars. Half arc typically between 450 
and 725 pc. This prediction may rise to 10-15 stars if one includes an age-dependence to the 
X-ray luminosities (Preibisch & Feigelson 2005b) into the simulations, since ~ 1/4 of the 
stars have ages < 1 Gyr. 

The 13 Chandra sources highlighted in Figure 15 have NIR magnitudes and colors 
lying within the locus of expected foreground stars. The spatial distribution of the 13 soft- 
spectrum stars is relatively uniform in the ACIS field (see the red x symbols in Figure 13b) 
and they have soft X-ray spectra expected from MS coronal emission (sources with flag 'f ' 
for 'foreground' in Table 3). We omit these X-ray sources from further consideration of the 
Cepheus population. 



5.3. Galactic background contamination 

The Besangon model is again convolved with X-ray luminosity functions to simulate 
the number of stars behind the cloud that may enter our X-ray sample. We consider only 
plausible classes of X-ray emitters (e.g., omitting MS A stars and subgiants) and adopt a 
limiting magnitude V < 22. The model then predicts ~ 1700 F dwarfs, ~ 550 G dwarfs, 
~ 1950 K dwarfs, ~ 25 M dwarfs, and ~ 632 giants. We again use the dwarf XLFs estab- 
lished for the solar neighborhood (Schmitt et al. 1995; Schmitt 1997). For giants we assume 
the lognormal XLF of ~ <28.9> ±0.9, obtained from Table 2 of Pizzolato et al. (2000) 
for a sample of 120 post-main-sequence late-type stars within 100 pc. Rarely, simulations 
result in up to 5(7) dwarfs (giants) detected in our Chandra observation, while typical runs 
suggest 0(3) dwarfs (giants). As we can not make individual determinations for 3 possible 
background giants among our PMS stars, we ignore Galactic background contamination in 
our further consideration of the Cepheus PMS population. The simulations yield the total 
typical observed X-ray fiux of the Galactic background population of up to Fx ~ 2 x 10"^"' 
erg cm~^ s~^ that may be contributing to the X-ray diffuse emission seen in the ACIS-I field. 
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6. Discovery of two rich stellair clusters 

Cephcus B, the hottest CO component of the Cepheus molecular cloud, is located at 
the northwestern edge of the cloud near the Cepheus 0B3 association. The association is 
thought to be composed of two subgroups of different ages, with the youngest lying closer 
to the molecular cloud (Sargent 1979). The interface between the molecular cloud and the 
OB association is delineated by the optically bright HII region S 155 with ionization fronts 
bounding the dust /molecular cloud (Felli et al. 1978). Testi et al. (1995) detected four radio 
structures — one extended source tracing an ionization front and three compact sources 
inside the molecular cloud — accompanied by three very red near-infrared objects. They 
suggested these are manifestations of a previously unrecognized young embedded stellar 
cluster at the interface of the HII region and the molecular cloud. 

Moreno-Corral et al. (1993) claimed the discovery of more than 50 NIR sources, with 
60% of them located close to or within the Cep B molecular cloud. Unfortunately, it seems 
that positions for those sources reported in their Table 5 are incorrect, as none of them 
matches known 2MASS or Chandra objects considering any possible systematic shift. Thus 
we omit the results of this work from further consideration. 

Naylor & Fabian (1999) found more than 50 ROSAT X-ray sources in Cep 0B3b; 
these revealed a low-mass stellar population accompanying the high-mass stars of the OB 
association. Optical CCD photometry and spectroscopy identified 10 ROSAT sources in 
their catalog, of which 4 were confirmed as T Tauri and 1 was suggested as a possible PMS 
star (Pozzo et al. 2003). Six (five) more optical stars were associated with T-Tauri (possible 
PMS stars) totaUing 10 (6) confirmed T-Tauri (possible PMS). The results of Testi et al. 
(1995), Naylor & Fabian (1999) and Pozzo et al. (2003) are summarized in Figure 16 which 
shows the previously known PMS stars lying in the ACIS-I field: 4 radio (boxes); 3 red 
NIR sources (x); 27 out of 50 ROSAT (marked as diamonds)^; and 5 out of 16 optical 
PMS stars (-I-). The four stars from Pozzo et al. (2003) with Chandra counterparts are 7a 
(Chandra #193), 8 (# 156), 12 (# 197), and 106 (# 160). Pozzo et al. (2003) recognized the 
incompleteness and low spatial resolution of the ROSAT data and anticipated our Chandra 
X-ray observation which would be capable of reaching a much larger PMS population. 

The popular scenario for star birth in the S 155/Cepheus B region (Minchin et al. 1992) 
is that the ultraviolet (UV) radiation from the 07n star HD 217086, the brightest member 
of Cep 0B3b chister, ionizes ambient material and induces a burst of star formation. We 



■^AU except one of these ROSAT sources has one or multiple Chandra counterparts. ROSAT source # 
41 is missed, probably because it lies on a gap between CCD chips in the ACIS array. One of the candidate 
optical PMS stars similarly lies on a chip gap. 



- 15 - 



investigate this model by searching for the next (third) generation of young stars directly 
behind the ionization front. Figure 13c shows the A^j^-stratified distribution of the Chandra 
source population, excluding field and extragalactic candidates, overlayed on an outline of 
the optical HII region. 

A clear congregation of ~ 64 highly-absorbed YSOs is seen in a ~ 3' (~ 0.6 pc) radius 

area, delineated by the black circle, at the Hll-molecular cloud interface. The remaining ~ 
321 Chandra stars with moderate or low absorption are located outside the molecular cloud, 
although ~ 9 highly absorbed stars are dispersed throughout the field. Figure 13d shows 
the distribution of Chandra stars with a JT-band excess inferred from 2MASS photometry 
(Figure 4). These stars are likely to be younger than the stars without bright inner dusty 
disks. These are preferentially associated with the cluster of highly- absorbed Chandra stars: 
26% of stars within the black circle have X-band excesses compared to 4% elsewhere in 
the field. We recognize that the UV radiation from the 07n star could evaporate nearby 
circumstellar disks, but this likely affects only the outer cooler disk regions. 

Our findings support the triggered star formation scenario of Minchin et al. (1992) with 
an unambiguous discovery of the third-generation embedded star cluster forming close to 
the ionization front. This is based on the spatial concentration of highly- absorbed Chandra 
sources and their high proportion of i^-band excesses just behind the S 155 ionization front. 
It seems unlikely that this cluster was formed by spontaneous self-gravitational collapse of 
molecular cloud material. Wc detect ~ 64 members of this embedded chistcr (henceforth 
denoted as "Ccp B embedded cluster"), which arc marked "e" (= "embedded") in column 
13 of Table 3. The 321 members of the older second generation cluster (hence denoted as 
"Cep 0B3b cluster within ACIS-I field" = "Cep 0B3b unobscured cluster") are marked "1" 
(= "lightly absorbed"). 



7. XLFs, KLFs and IMFs 

A major goal of studies of young stellar clusters is to quantify and understand the 
distribution of masses arising from star formation processes (Corbelli et al. 2005). This 
effort is subject to a variety of difficulties, but as X-ray surveys present different problems 
than those arising at other wavebands they can be valuable for IMF studies (Feigelson & 
Getman 2005). In regions like Cepheus B, NIR surveys are overwhelmed by unrelated field 
stars (e.g. Figure 4). A common procedure is to subtract statistically a NIR source counts 
from a nearby control field and convert the resulting i^T-band luminosity function (KLF) into 
masses using theoretical stellar models to estimate a cluster IMF (Lada & Lada 1995). We 
use a similar procedure on the Cepheus B region below (§7.3). X-ray surveys, on the other 
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hand, suffer relatively little contamination; we found that no more than 5% of Chandra 
sources in the Cepheus B field are likely extragalactic contaminants and no more than 4% 
are likely stellar contaminants (§5). Infrared surveys also are biased towards selection of 
young stars with heavy protoplanetary disks, while X-ray surveys arc mostly unbiased with 
respect to disks (although a small bias against accretion is present, Preibisch et al. 2005). 
X-ray surveys penetrate deeper into obscuring material than optical or NIR surveys, but are 
less complete in detecting the lowest mass objects. 

In this section, wc first examine the empirical stellar XLF of the Cepheus population 
and compare it to the well-studied Orion Nebula Cluster (ONC) XLF (§7.1). Wc then use 
independently derived masses for Cepheus members to estimate the IMF (§7.2) and compare 
this with an IMF created from the KLF (§7.3). The result is surprising: the Cep 0B3 cluster 
— at least the part of it seen in the ACIS-I field — appears to have a deficiency in higher 
mass stars compared to the ONC and standard IMFs. 

7.1. Cep OB3b and ONC XLFs 

Comparison of XLFs between poorly characterized young stellar clusters like Cep 0B3b 
and a well-characterized cluster like the ONC can give insight into any missing population 
of stars. This was done by Grosso et al. (2005) to estimate the total population of the 
embedded Beckhn-Neugebauer-Kleinman-Low OMC-1 population. Figure 17 shows such 
a comparison among the unobscured COUP population (iipper solid line), the Cep 0B3b 
population (dashed line with error bars), and the Cep B embedded population (bottom solid 
line with error bars). The ONC sample (= 839 COUP lightly obscured cool stars) consists 
of COUP cloud members excluding OB stars and sources with absorbing column density 
log Nh > 22.0 cm~^ {Ay > 5 mag) (Feigelson et al. 2005). XLF error bars for the Cep 0B3b 
and Cep B embedded populations show the range containing 68% (1 a equivalent) from 
Monte Carlo simulated XLF distributions where individual X-ray luminosities are randomly 
drawn from Gaussian distributions with mean equal to the measured source's logLx and 
variances based on observed AlogLx values. The error bars thus include both the VN 
counting error and the Lx-dcpcndent measurement error. Regression lines to COUP's XLFs 
were added to figures by eye (red solid lines) and scaled vertically to obtain a predicted XLF 
for the Cep 0B3b population (red dashed lines). 

It is clear from Figure 17 that the Cep 0B3b cluster XLF shows a significant excess of 
stars around logL^ c ~ 29.7, log Lt ~ 29.5, and log Lh^c ~ log Lh ~ 29.3 compared to the 
scaled ONC XLF. This excess is not related to incompleteness of the current observation. 
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Figure 18^° shows that the completeness hmit for the observed total band luminosity is 
around logLj = 29.5 erg s~^, and is still above 50% at logLj = 29.0 erg s~^. This is far 
below the range logL^ > 30.2 erg s~^ over which the ONC and Cep 0B3b XLFs were scaled. 
We arc thus confident that the relative amplitudes of the XLFs in Figure 17 arc correct. In 
fact, the excess of Cep 0B3b X-ray stars over the ONC XLF becomes even stronger if the 
completeness curve from Figure 18 is applied to the incomplete Cep 0B3b Lj-XLF of Figure 
17. 

Table 5 summarizes total X-ray luminosities for all X-ray detected stars in the Cep 0B3b 
unobscurcd and Cep B embedded populations within the AClS-1 field. It thus appears that 
the Cep 0B3b population within the ACIS-I field has an XLF with a difi^crcnt shape than 
seen in the ONC. Instead of a two-powerlaw model with break around Lt^^ ^ 2.5 x 10^° 
erg s~^, the Cep OBSb XLF with the sensitivity limit of logLt^cUm ~ 10^^ erg s"^ can be 
roughly modelled as a lognormal with mean around < log L^^c >~ 29.8 erg s~^ and standard 
deviation of a(\ogLt c) ~ 0.6. The mean and standard deviation of cluster XLFs are perhaps 
observation biased - for example for more sensitive observations of NGC 1333 and IC 348 
with log Lt^cHm ~ 10^^ erg s~^ their XLF shapes are wider with the standard deviation of 
crilogLt c) ~ 0.9, peaking at lower luminosities of < log L^^c 29.5 erg s~^ (Feigelson & 
Getman 2005). 



7.2. Cep OBSb and ONC IMFs 

To compare the Cepheus IMF with that of the ONC, we need to adopt compatible 
procedures for estimating masses. To be consistent with the current work and to derive 
masses for the whole sample, we re-derive^^ COUP stellar masses here using NIR photometry 
(Table 10 in Getman et al. (2005)). Similar to Cepheus masses, ONC mass estimates are 
based on star locations in the J vs. J — H color-magnitude diagram and theoretical stellar 
model tracks derived by Baraffe et al. (1998) (for 0.02 Mq ^ M ^ 1.4 Mq) and Siess et al. 
(2000) (for 1.4 Mq ^ M ^ 7.0 Mq). We are aware that these mass estimates are subject to 



^^Figure 18 is derived from Monte-Carlo simulations performed by randomly drawing sources of a given 
luminosity scattered throughout the ACIS image and applying local background levels found in the image, 
similar to the AGN simulations in §5.1. Observed total band luminosity Lt was chosen as it is more readily 
converted to the observed source count rate assuming a typical distribution of source median energy, and 
the source significance then is easily calculated. Source significance thresholds are then evaluated to mimic 
the behavior of our source detection procedure. 

Table 9 in Getman et al. (2005) provides masses for half of the COUP sample based on optical spec- 
troscopy. 
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significant uncertainties and may be incompatible with masses obtained by other methods 
such as optical spectroscopy. However, here we are not so much interested in exact masses 
for individual objects as in the general trends of mass distributions and consistent treatment 
for comparison between clusters. 

The IMF for the COUP 839 sample is shown as an upper dashed histogram in Figure 
19. The COUP sample can be extended to account for incompleteness from instrumental 
causes: low local ACIS exposure time, high local background due to readout streaks from 
bright ONC sources, and contamination from the point spread function wings of close bright 
X-ray sources. There arc 82 such ONC stars with Ay ^ 5 mag and Ks < 13 mag, corre- 
sponding to M > Q.l Mq in Orion (Getman et al. 2005b). The resulting IMF of the 921 
(839-1-82) extended COUP ONC sample is shown as an upper solid histogram on Figure 19. 
Comparison with the ONC IMF fit based on KLFs obtained by Muench et al. (2002) (dashed- 
dotted hne) indicates that the sample is complete down to 0.1 Mq and about 50% complete 
(120 versus ~ 240 stars) within the low-mass stellar-sub-stellar range of (0.03 — 0.1) Mq. 
This incompleteness at low masses is an obvious consequence of an interaction between the 
observational flux limit and the X-ray/mass correlation (Preibisch et al. 2005). 

We now turn to the Cep 0B3b IMF derived in a similar manner. From 321 (64) PMS 

stars of the Chandra Ccphcus unobscured (obscured) population, 310 (55) have masses 
estimated from NIR photometry (Table 2). Figure 20 confirms that the PMS stars in Cepheus 
(red - unobscured, blue - obscured) follow Lx versus mass correlations very similar to that 
seen in COUP stars (grey dots). Even observed luminosity Lt, prone to changes as absorption 
varies, has a similar correlation with mass, as the median column densities for COUP and 
Cepheus unobscured populations are close to each other. Constant and slope parameters 
derived from the EM linear regression fit for the COUP unabsorbed sample are given for 
all four measures of X-ray luminosities in Table 6. The regression parameters are slightly 
different from those derived in Preibisch et al. (2005) for the ONC lightly absorbed optical 
sample, but this can arise from the higher sensitivity of the COUP 839 sample than that 
of the optical sample, and from different methods used to estimate masses (Preibisch et 
al. 2005). For 11(9) stars of the Cepheus unobscured (obscured) populations with unknown 
masses, we can use their measured X-ray luminosities with the regression fit parameters from 
Table 6 to get rough estimates on their masses and to add them to the appropriate bins of 
the IMF histograms. This is done in the lower histograms of Figure 19. 

Though the histogram for the Cep 0B3b population within the ACIS-I field (thick 
dashed line) is rough and clearly incomplete below ~ 0.2 Mq, we see a good agreement with 
the ONC IMF from 0.5 to 3 Mq but an excess of ~ 20 or more sources (depending on the 
choice of scahng the ONC-Muench IMF) in the log(M/Mo) 0.5 (0.3 Mq) bin. Expressed 
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in another way, the Cep 0B3b IMF in contrast to the ONC does not exhibit the turnover 
from the Salpeter IMF around 0.5 - 0.6 M©. 

7.3. Cep OB3b KLFs 

We can estimate the total population of the Cepheus stellar clusters by extrapolating 
the IMFs below the incompleteness limit produced by the Chandra flux limit. 

We obtain a coarse statistical sample of 2MASS PMS members of the Cep 0B3b popula- 
tion in the ACIS-I field by subtracting a local background of foreground and background NIR 
objects from a control field of the 2MASS catalog. An IMF complete down to M = O.OSM© 
at a distance of 725 pc corresponds to J ~ 16 mag which is roughly the 2MASS completeness 
hmit. Figure 13(c) suggests that ~ 90% of the Cep 0B3b population are clustered within 
an area of ~ 127 arcmin^ lying within the ACIS-I field. This area roughly matches the 
"pentagon" shown in Figure 21 where X-ray extended emission is seen (§8). We extracted 
2MASS sources from 30 control fields with 127 arcmin^ area lying in a 3° x 3° neighborhood 
of the Cep B cloud as shown in Figure 22. It is not trivial to establish that the control fields 
are free from PMS stars, as the entire region is occupied by the members of the sparse stellar 
0B3 association and we trace a gradient in NIR source density as one progresses away from 
the cluster core. Towards the northwest extracted 2MASS source density clearly reflects the 
evolution of PMS source distribution in the Cep 0B3 association — the number of sources 
systematically decreases down to 830 objects (at about 45' (~ 9.5 pc) away from the ACIS-I 
field), but then systematically rises again. The area with the least number of objects, 830 in 
a 127 arcmin^ circle centered on {a, 6) = (342. 92539, +63. 26147), perhaps lies at the bound- 
ary between different sub-groups of the Cep 0B3 association ("N-W" background area); see 
Figure 1 in Sargent (1977). A second similar area is found toward the west, 814 2MASS 
sources at (341.68270,+62.27262) as shown in Figure 22 ("W" background area). Figure 23 
compares cumulative distributions of J — if color excesses of 2MASS sources seen towards 
the two control fields (dashed fines) with that seen towards the molecular cloud (dotted line) . 
The large population of (J — if) < 1 sources in the control fields supports the idea that they 
are relatively unobscured areas unaffected by the Cepheus molecular cloud. 

To obtain a coarse IMF trend of the 2MASS unobscured sample of PMS stars within 
the ACIS-I field we use the correlation Kg — 10.86 — 3.32 x log(M) based on the theoretical 
synthesized 1 Myr PMS model, Siess -|- Baraffe, employed here for < Ay >~ 2.5 mag 
{Ak ~ 0.28 mag). This permits conversion of the KLF to stellar mass and production of 
two mass histograms (black solid lines in Figure 24a) after subtraction of the two control 
fields with accounting for 10% of the unobscured Cep 0B3b members dispersed outside of 
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the "pentagon" area. 

The result is that our Chandra unobscured sample of 321 stars (thick dashed in Figure 
24a) appears complete down to ~ 0.4 Mq, and that 400 — 450 stars are missing from the 
Chandra sample between 0.08 and 0.4 Mq. Figure 24a also compares this estimated IMF 
with that obtained for the ONC by Muench et al. (2002). 

7.4. Discussion 

Comparison of the IMF based on the 321 X-ray detected PMS stars with the IMF based 
on the background-subtracted KLF suggests that the Chandra image is complete down to 
~ 0.4 M0 (Figure 24a). The difference between the infrared and X-ray source counts implies 
that 400 — 450 cluster stars are missing from the Chandra image, giving a total unobscured 
cluster population of ~ 720 — 770 stars in the ~ 130 — 170 arcmin^ area. We estimate the total 
X-ray luminosity of these undetected stars by multiplying the number of undetected stars 
by a typical X-ray luminosity for the appropriate mass bin using regression formulae from 
Table 6. Table 7 summarizes results of this exercise and shows that a total of logLj ~ 31.3 
and log Lt^c ~ 31.5 ergs s~^ is expected from the undetected population. 

Assuming that the IMF for the embedded cluster also follows the mass distribution of 
the unabsorbed cluster, the embedded cluster of 64 X-ray-detected members is complete 
down to M = 0.8 Mq. We then estimate that 300 — 340 undetected obscured stars down to 
M ~ 0.08 M0 are hidden in the cluster, giving a total of ~ 360 — 400 obscured stars. 

We find that both the XLF and IMF comparisons indicate that the Cep 0B3b population 
within the AClS-1 field exhibits, compared to the ONC, more stars around M ~ 0.3 M0 
relative to 0.5 < M < 3 Mq stars. This difference cannot be attributed to incompleteness in 
the COUP observation, as the addition of 82 stars missing from COUP does not significantly 
change the shape of COUP's IMF and hence XLF, and will not remove the discrepancy. 
It also cannot be attributed to incompleteness in the Cep 0B3b observation (Figure 18); 
inclusion of missing stars below log Lt = 29.0 erg s~^ would only increase the inferred number 
of lower mass cluster members. Similarly, the inclusion of binaries unresolved in Cep 0B3b 
but resolved in the ONC (Getman et al. 2005b) would enhance the higher, not the lower, 
ends of the Cep 0B3b XLF and IMF. 

It is possible that our assumption of similar {-^ 1 Myr) ages for all PMS stars of the 
Cep 0B3b population within the ACIS-I field is incorrect. Recall that Pozzo et al. (2003) 
found the age of ~ 1 Myr for four out of six and > 3 Myr for another two optical PMS 
candidates in our ACIS-I field. For the larger Cep 0B3b association outside of the ACIS-I 
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field, the observed age spread on the V vs. {V — I) color-magnitude diagram becomes even 
more dramatic ranging from < 1 Myr to 10 Myr, and is too broad to be accounted for by 
binarity, photometric uncertainty and variability (Burningham et al. 2005). Subgroups of 
different ages may exist, each with different XLF and IMF characteristics. The XLFs are 
relatively insensitive to these age effects; between 1 and 10 Myr, only a slight decay (roughly 
0.3 in log Lt) is found in the COUP population (Preibisch & Feigelson 2005b). An older 
assumed age will have a bigger effect on the KLF. Figure 24b summarizes a KLF analysis 
similar to that in §7.3 but assuming a cluster age of 3 Myr. In this case. X-ray detected 
stars are more massive with the completeness hmit shifted from ~ 0.4 M© (for 1 Myr) down 
to ~ 0.8 Mq (for 3 Myr). The excess of lower mass stars over the ONC IMF then mostly 
disappears. However, a younger age for Ccp 0B3b has two advantages: the unusual inferred 
IMF agrees with the unusual observed XLF (§7.1), which is insensitive to age assumptions; 
and the shape of the background-subtracted KLF more closely agrees with the IMF of the 
X-ray detected population for 1 Myr compared to 3 Myr (Figure 24). An argument in favor 
of the older age can be seen in Figure 24b: the background-subtracted KLF more closely 
matches the Galactic field IMF (grey line) than in the case of 1 Myr. 

We thus emerge with a clear determination that the XLFs of Cep 0B3b and the ONC 
differ, but without a clear interpretation of this effect. It may be due to a different IMF with 
more lower-mass stars in Cep 0B3b, or to differences in age and history of star formation 
in the two clusters. FoUowup optical photometry and spectroscopy of 321 Cep 0B3b X-ray 
sources would assist in unravelling these possibilities, as individual (rather than statistical) 
masses and ages can be directly estimated on the Hertzsprung- Russell diagram. 

If we assume that the 1 Myr ages are correct, then the inferred IMF difference may 
be related to differences between the physical processes of star formation in the Cep B and 
Orion molecular clouds. The average Jeans masses may be lower in Cep B than in Orion due 
to different densities or pressures. However if the Jeans mass is universal around Mj ~ 0.1 
M0 (Elmegreen 1999), then the IMF difference could arise the prevalence of different star for- 
mation mechanisms in the two clouds (Elmegreen 2004). In this multi-component scenario, 
lower-mass stars would have formed similarly in Cep B and Orion from pre-stellar conden- 
sations of Mj mass, but higher-mass stars would have grown from coalescence of pre-stellar 
condensations and gas accretion more efficiently in the more crowded Orion environment. 
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8. The diffuse X-ray emission 

Soft X-ray diffuse emission is clearly present in the central part of the AClS-1 field. 
Figure 2f shows an adaptively smoothed map of the field with all 43f point sources removed^^. 
The diffuse emission occupies the pentagonal f 27 arcmin^ area where 90% of the clustered X- 
ray detected members of the Cep 0B3b population reside (§7.3). We extracted the spectrum 
of the diffuse emission from the pentagonal area (see footnote 12) and subtracted extraneous 
spectral components associated with the instrumental, part of cosmic X-ray, and field star 
backgrounds measured on the molecular cloud in the southeastern corner of the field, as 
shown in Figure 21. The resulting spectrum of the diffuse emission is shown in Figure 25(a). 
It can be reasonably well fit by a single-temperature optically thin thermal plasma with 
kT ~ 0.8 keV subject to logiY^ ~ 21.7 cm'^ (Ay ~ 2.5 mag) absorption, which is typical 
for the unobscured Cepheus region. The observed (absorption-corrected) X-ray luminosity 
of the diffuse emission is logLj ~ 31.3 (logL^^c — 31.8) erg s~^. 

This is substantially weaker than the log Lf ~ 32 — 33 erg in diffuse X-rays attributed 
to colliding O star winds in other HIT regions such as M 17 (Townsley et al. 2003). In theory 
a strong wind from a single massive star, interacting with the surrounding ISM, may fill 
the surrounding volume with hot shocked gas and produce an X-ray emission as high as 
1035-36 gj,g g-i rj.^^ 07n HD 217086 star has a high wind terminal velocity of Ko = 2510 
km/s derived from ultraviolet spectroscopy by Prinja et al. (1990). Such a fast wind may 
heat the gas to the temperatures of 10^~*^ K. On the other hand, HD 217086 has a low 
mass-loss rate of M ~ 2 x 10~^ Mq yr~^ (Howarth & Prinja 1989). Stevens & Hartwell 
(2003) have calculated a scaling between expected X-ray emission and the wind scaling 
factor Xci = M'^/{RcVoo) at temperatures higher than 10^ K. For a cluster core radius of 
i?c ~ 1 pc and the HD 217086 wind parameters, the predicted diffuse X-ray luminosity is 
logLx < 28 erg s~^ for the wide range 0.1 — 1 of the thermalization efficiency parameter 
T). This is at least 3 orders of magnitude smaller than the X-ray luminosity derived for the 
detected diffuse X-ray emission (Table 8). This casts considerable doubt on the idea that 
the O star wind shocks are responsible for the emission. 

The alternative and more natural scenario for the origin of the diffuse emission in 
Cepheus is X-ray emission from the 400 — 450 lower mass PMS stars, which are not individu- 
ally detected (§7.3) due to the high sensitivity hmits {logLt^cUm — 29 for the hghtly-obscured 



^^Details on the production of such images, and extraction source-diffuse and background spectra can 

be found on-line in the recipes to ^i? http://www.astro.psu.edu/xray/docs/TARA/ac_uscrs_guidc/, §7.1.2. 
They have been previously applied to Chandra images of high mass star formation regions by Townsley et 
al. (2003). Appendix C of Townsley et al. (2003) describes the method. 
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population; Figure 17) of the current Chandra observation. Recall that the Galactic back- 
ground contribution to the diffuse emission is small, log L-t < 30, as indicated by the simula- 
tions described in §5.3. Several arguments support the explanation for the diffuse emission 
from an unresolved PMS component in Cepheus. 

1. The diffuse emission spatially coincides with the clustered X-ray detected members of 
the lightly-obscured Cep 0B3b population as expected if cluster stars are responsible. 
Compare the distribution of detected stars in Figure 13 with the diffuse emission in 
Figure 21. 

2. The spectral shape of the diffuse emission in Cepheus closely resembles the composite 
spectrum of the detected X-ray population as shown in Figure 25 for the Cepheus ob- 
servation. A similar finding can be seen in the cool lightly obscured COUP population 
(839 sources introduced in §7.1). Figure 26 illustrates the resemblance between (panel 
a) a COUP sub-sample containing the weakest (log Lt^c < 29) sources — corresponding 
to the Cepheus diffuse emission here attributed to the undetected population — and 
(panel b) the full COUP sample (820 sources, excluding those piled-up) — correspond- 
ing to the detected Cepheus population. The plasma temperature of kT ~ 0.8 keV is 
similar for both the Cepheus diffuse emission and the COUP weak sources (Table 8). 
This is not a surprise, as it has been established from the COUP project that the tem- 
perature of the hot plasma component for ONC T-Tauri stars having two-temperature 
spectral fits decreases as X-ray luminosity drops, while the cooler plasma component 
remains around 10 MK (~ 0.86 keV) in nearly all T-Tauri stars (Preibisch et al. 2005). 

3. The inferred integrated (absorption-corrected) X-ray luminosity of the missing PMS 
population of the unobscured Cep 0B3b cluster, estimated in §7 to be logL^ ~ 31.3 
(logLf c — 31.5) erg s^\ is very close to the observed diffuse X-ray luminosity of 
log Lt ~ 31.3 (logLi^c — 31.8) erg s"^.-*-^ 



^^The difference of a factor of ^ 2 between absorption-corrected X-ray luminosities of log Lt^c — 31.5 
expected from faint cluster stars and 31.8 erg s~^ observed in diffuse emission can be explained as a statistical 
artifact. The same discrepancy is seen for the integrated absorption-corrected luminosities of 326 COUP weak 
sources considered individually in comparison with the absorption-corrected luminosity derived from their 
composite spectrum. When considered individually, the poor statistics of these sources leads to successful 
fits by one-temperature with median plasma energy kT ~ 1.6 keV. This is noticeably harder than the 
kT ~ 0.8 kcv energy derived from the combined fit. This arises because of a bias towards detecting photons 
at energies where the Chandra-AClS effective area is greatest around 1.5 keV; substantially softer or harder 
photons are not detected in sufficient quantity in individual weak sources to be represented in the XSPEC 
fitting procedure. Due to this bias, absorption-corrected X-ray luminosities Lt,c of weak sources may often 
be underestimated due to the unseen soft temperature component, as discussed by Getman et al. (2005); 
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A strong Ne IX and Ne X line complex is present in the composite spectrum of the 
Cep 0B3b unobscured PMS members (Figure 25(b)). This emission line is prominent in the 
composite spectrum of COUP ONC stars (Figure 26(b)), and Feigelson et al. (2005) have 
argued that this line should be a useful diagnostic of magnetically active PMS stars. 



9. Two interesting X-ray stars 

9.1. Wind shocks in HD 217086 (07n) 

Among the members of the Cep 0B3b population within the AClS-1 field, only one 
massive star is seen in X-rays: the 07n star HD 217086, Chandra source #240. It is 
thought to be the principal ionizing source of the S 155 HII region. 

Our understanding of the production of X-rays in the inner winds of O stars is somewhat 
confused. Some stars show the soft, constant X-ray emission expected from a myriad small 
shocks in unstable radiatively accelerated winds (Lucy & White 1980). But others show hard 
spectral components, rotational modulations, flaring light curves, and unexpected spectral 
line shapes that indicate other processes are operative. Stelzer et al. (2005) performed an 
in-depth study of the X-ray emission of the early-type stars of the ONC and found that 
only two of nine hot massive stars with spectral class 07-B3 exhibit exclusively the constant 
soft-spectrum at \og{Lx/Lhoi) ~ ~7 expected from the standard model involving small- 
scale shocks from unstable radiatively accelerated winds. Most of the other massive stars, 
including the brightest 9^ Ori C (07), show combinations of soft wind emission and other 
phenomena which could be explained by a magnetically confined wind shock model, where 
the harder and variable X-rays are produced by magnetically-mediated large-scale shocks 
(Babel & Montmerle 1997; ud-Doula & Owocki 2002). 

Though the Chandra source ^^240 falls by chance on an ACIS CCD chip gap, this should 
have little effects on its properties once we evaluate the exposure map value at source's 
location. Its effective exposure is 12.9 ks (Table 1), about 2.3 times shorter than the 30.1 ks 
effective exposure for the whole observation. The source spectrum and lightcurve are shown 
in Figure 27. To ensure that the interaction of sateUite dithering with the chip gap does not 



Feigelson et al. (2005). 

^^One may note strong residuals between the spectral model and data around 2.1 keV in (Figure 25(b)) and 
26(b). These are likely instrumental features, related to the fact the effective area of the Chandra telescope 
mirrors is not properly calibrated at the iridium M edge near 2.1 keV, possibly due to a thin contaminant layer 
on the mirrors. For more details see http://cxc.harvard.edu/ccw/proceedings/04_proc/presentations/jerius/. 
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affect the shape of the hghtcurve, the hghtcurve is shown with 2 ks binsize, which is two times 
larger than the dithering period. A Kolmogorov-Smirnov test shows only marginal evidence 
{Pks ~ 0.07) for variability during the 8 hour observation; it was similarly constant during 
the earlier ROSAT observation (Naylor & Fabian 1999). The X-ray spectral analysis gives 
the count rate of ~ 32.3 cnt ks^^, plasma temperature ~ 0.6 keV, and the X-ray luminosity 
in (0.5-8.0) keV band of log L^^c ^ 31.5 erg s~^ after correction for log Nh — 21.7 cm~^ 
absorption. Comparison with the coarse bolometric luminosity derived from the 2MASS 
color-magnitude diagram (Table 2) gives log{Lt^c/ Lboi) — —7 with a stellar mass of ~ 24Mq. 

Thus its soft constant X-ray emission at \og{Lt,c/ Li,,,}) = —7 level is probably at- 
tributable to the traditional Lucy- White model of small-scale wind shocks without invoking 
magnetic fields. 

9.2. Superflare in Chandra star # 294 

The brightest X-ray source in the ACIS field is source ^^294 with 2938 counts extracted 
from the 89% encircled energy region. Source # 294 has an extraordinarily luminous and 
hard X-ray flare whose amplitude, power, and temperature are similar to the brightest flares 
seen in the Orion population (Grosso et al. 2004; Favata et al. 2005). Its hghtcurve has the 
typical morphology of stellar flares with the fast rise and the gradual decay (Figure 28(a)). 
During the ~ 1.1 hour rise phase the count rate rose by a factor ~ 20, ~ 50, and ~ 10 in the 
full, hard, and soft bands, respectively. Due to the short exposure of our observation, only 
the flrst ~ 4.4 hours of the decay phase were observed. We divide this phase into a rapid 
and a slow cooling phases. 

Spectra extracted from the quiescent, rise and peak phases, rapid cooling, and slow 
cooling phases are compared in Figure 28(b). They all can be successfully represented 
by single-temperature optically thin plasma models seen through a similar amount of gas 
absorption of log Nh = 21.8 cm~^. During the quiescent phase, the plasma had an average 
temperature of ~ 1.7 keV (~ 20 MK) which heated rapidly during the flare rise phase of ~ 5 
ks up to <kT> > 15 keV. The nominal XSPEC flt gives 64 keV (- 740 MK). Over the 
next ~ 4.4 ks, the plasma rapidly cooled down to <kT> ~ 12 keV (~ 140 MK), and cooling 
continued during the slow cooling phase to <kT> ~ 4.5 keV (~ 50 MK) over the last 10 ks. 
X-ray luminosities rose from a quiescent level of Lt^c — 9.7 x lO^'^ {Lh,c — 3.5 x lO^'^) ergs 

in the full (hard) energy bands to Lt^c = 2.2 x 10^^ {Lh^c = 1-8 x 10^^) erg s~^ during 
the rise-peak phase. These temperatures and luminosities are similar to the most luminous 
X-ray flares ever recorded in young stellar objects: V773 Tau in the Taurus-Auriga clouds 
with Lpk > 10^3 erg s'^ and Tpk > 100 MK (Tsuboi et al. 1998), LkHa 312 in the M 78 
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reflection nebula of Orion with Lpk ~ 10^^ erg s ^ and Tpk ~ 90 MK (Grosso et al. 2004), 
and COUP #1568 in the Orion Nebula Cluster with Lpk ~ 7.8 x lO^^ gj-g g-i j.^^ ^ 599 
MK (Favata et al. 2005). 

With detailed data on the decay of flux and temperature during the cooling phase, we 
can infer properties of the radiating plasma structures within the framework of the time- 
dependent hydrodynamic model in a coronal magnetic loop. Reale et al. (1997) establish 
a formula for estimating loop size with a proper account for the possibility of prolonged 
heating during the decay phase. They find L oc Tic^/Tpk/ F{Q, where L is the length of the 
loop, Tic is the flare decay time, and Tpk is the plasma temperature at the peak. -F(C) is 
a correction factor for prolonged heating that is a function of the slope C, of the trajectory 
in the temperature-density diagram. F{Q must be calibrated for each X-ray telescope and 
detector. The model assumes the plasma is confined within the semicircular loop of uniform 
cross-section, whose geometry remains unaltered during the whole flare, and the motion and 
energy transport occur along magnetic field lines of the loop. Favata et al. (2005) applied 
this approach to infer the size of the fiaring structures for bright X-ray flares in the Orion 
Nebula cluster population with F{Q calibrated to Chandra and ACIS-I. They found several 
superflares with inferred loop sizes long enough to link the stellar photosphere with the inner 
rim of the circumstellar disk. 

The slope C, is usually measured in the logT vs. \og \J\EM) plane, where ^J\EM) is 
used as a proxy for the plasma density. In the case of the source # 294 flare, we adopt 
a slightly different technical approach. Instead of performing time-sliced spectroscopy with 
XSPEC over a few characteristic time blocks of the lightcurve as shown in Figure 28(b), we 
employ an adaptivcly smoothed median energy time series to derive hundreds of time points 
of {kT, EM) along the decay phase of the hghtcurve. This is achieved by calibrating median 
energies to temperatures using the spectral simulations described in §4. Median energies 
are more statistically stable for small- datasets than kT parameters of nonhnear XSPEC 
regressions. We assume Nh is constant and we ignore the weak soft spectral component that 
is likely present from the quiescent background. 

We test this method on the decay phase of the flare in the Orion COUP source #1343 

studied in detail by Favata et al. (2005). The lightcurve for the COUP source #1343 is shown 
in Figure 29 with black dots representing the lightcurve histogram, binned with a bin size of 
~ 1 hour. The peak-decay phase studied in our analysis ranges from 230 to 330 ks from the 
beginning of the COUP observations. The black and grey lines represent smoothed versions 
of the lightcurve using the adaptively smoothed algorithm from AE with different minimum 
number of counts in the kernel (625 and 100 count minimum, respectively). The 100 count 
kernel is too small as it leaves many noise spikes, so we use the 625 count kernel although it 
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has oversmoothed the narrow peak phase of the hghtcurve. Results of the analysis appear 
in Figure 30. The top left panel shows the trend of the median energy with 88 time points 
over the 100 ks time range. The top right and bottom left panels show the derived decays of 
temperature and emission measure assuming a constant absorption of logNn = 21.9 cm~^. 
Finally the bottom right panel gives the logT - log ^/{EM) dia gram over the 100 ks peak- 
decay phase; the flare evolves from the upper right to the lower left. The best fit decay 
slope is C = 1-86 ± 0.29. These last three panels look very similar to those shown in Figure 
3 of Favata et al. (2005) who obtain a similar slope of C = 1-95 ± 0.51 using coarser time 
samphng and XSPEC fits. This validates our use of adaptively smoothed median energies 
as estimates of plasma temperatures in flare loop models. 

The method is apphed to 12.3 ks (15.5 — 27.8 ks after the start of the observation) of 
the ACIS-I source # 294 lightcurve; results are shown in Figure 31. Trends of the median 
energy for two different smoothing kernels are alike within the studied decay phase of the 
flare. We report results for the 570 count kernel, although it does a poor job tracing the 
earlier rapidly changing flare peak. The derived temperatures are very high, at least T ~ 300 
MK and conceivably as high as 800 — 1000 MK during the first 1 ks of the decay phase. But 
they rapidly drop to T ~ 100 MK over the next hour and then decrease more gradually to 
T" ~ 60 MK over the following 2 hours. 

The derived log T — log ^J\EM) correlation (bottom right panel) suggests two decay 
phases: rapid cooling over the first ~ 4.4 ks with = 8.08 ± 0.96, followed by slow cooling 
with C = 1.4 ± 0.13. The median temperatures in the upper right panels arc in good 
agreement with the cruder time-sliced spectroscopy shown in Figure 28(b). The slow decay 
value C = 1.4 is close to the upper limit of 1.5 for a freely decaying loop with no prolonged 
heating input (Favata et al. 2005; Reale et al. 1997). 

Since wc observe only the first ~ 4.4 hours of the decay phase, we can obtain only a 
lower limit to the phase duration. Tie > 15 ks, and thus only a lower limit to the inferred loop 
length. Using formulae (5)-(7) in Favata et al. (2005) and truncating our peak temperature 
around 300 MK, we derive a loop length L > 6.4 x 10^^ cm. We roughly estimate the stellar 
radius for source # 294 to be i?* ~ 2.5Rq using (Siess et al. 2000) PMS models for a ~ 1 Myr 
old star of ~ IAMq mass (Table 2). The magnetic loop is then > 3.7 times the radius of the 
star. This limit is even more conservative if the peak temperature exceeds 300 MK; the loop 
would be 2 times larger if a peak temperature around 800-1000 MK is used, as suggested 
from the plasma temperatures derived in Figure 31. Source #294 may thus be another case 
where magnetic fields connect the star to the circumstellar disk, as discussed by Favata et 
al. (2005). Though there is no an indication for the presence of K-band excess in source # 
294, it is well known that JHK photometry alone is not sufficient to detect circumstellar 
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material. Ground-based JHKLM or Spitzer Space Telescope photometry would allow much 
better discrimination between stars with and without circumstellar disks. 



10. Conclusions 

We present a 30 ks (net exposure) Chandra ACIS observation of the Cepheus B molec- 
ular cloud and a portion of the nearby Cep 0B3b OB association. The main conclusions of 
our study are as follows: 

1. We detect 431 X-ray sources with a limiting luminosity for the region outside of the 
molecular cloud of ~ 10^^ ergs/s. Three hundred eighty four of them are unambiguously 
identified with 2MASS NIR sources. Upon visual inspection of 2MASS atlases we find NIR 
counterparts for at least 10 (2) more soft (hard) X-ray sources, which have not been registered 
in the 2MASS catalog. On the NIR color-magnitude diagram most of the detected X-ray 
sources occupy the locus of low-mass PMS stars with a typical absorption of Ay ~ 2.5 mag 
(assuming age of ~ 1 Myr), but as many as 60 sources with NIR counterparts, residing close 
to/in the molecular cloud, have even much higher absorptions up to Ay ~ 28 mag. 

2. Twenty four unidentified hard X-ray sources are classified as possible cxtragalactic 
contaminants (Tabic 3, 'a'). Though some of them still may be embedded members of 
the cloud, about 70% percent of them are probably true AGNs, as argued by a detailed 
simulation of the cxtragalactic background population, including instrumental background 
and local (Cepheus) plus large scale Galactic absorptions. Through careful simulations of X- 
ray emission from the Galactic foreground and background populations, employing a stellar 
population synthesis model of the Galactic disk, and through the inspection of NIR colors 
of X-ray detected sources we expect to have no more than a few Galactic background stars 
and we suggest 13 field star candidates (Table 3, 'f'). 

3. In order to obtain a rehable X-ray luminosity function of the detected PMS stars 
we determine source X-ray fiuxes based on knowledge of their NIR absorptions and X-ray 
median energies. We then perform careful MARX simulations to derive errors on median 
energy, combine them with known errors on count rate and absorption, and propagate them 
to the flux errors through Monte-Carlo simulations. 

4. We discover more than 380 PMS stars in the ACIS-I field of the Cepheus region, 
which provide the best census to date for the stellar population of the cloud. Based on Nh 
and K-hand stratified spatial distributions of the detected PMS stars we separate them into 
two rich clusters: the lightly obscured cluster (= "Cep 0B3b within ACIS-I field") (321 stars. 
Table 3, '1') and the embedded cluster (64 stars. Table 3, 'e'). We estimate that another 
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400 — 450 (300 — 340) low mass PMS stars are present in the unobscured (obscured) clusters 
within the Chandra field but are too X-ray-faint to be individually identified. 

5. We find that the XLF of the hghtly-obscured Cep 0B3b cluster region within the 
ACIS-I field shows an excess of stars around logL^ ~ 29.7 erg s~^ (0.5 — 8 keV, after 
individual correction for absorption) compared to the well-studied Orion Nebula Cluster. 
This may be caused by a nonstandard IMF in Cep 0B3b with an excess of ~ 0.3 Mq stars 
if the cluster age is ~ 1 Myr. This excess also appears in the independently derived KLF 
based on 2MASS source counts (after an uncertain subtraction of background counts from 
two control fields) which show more ~ 13 — 14 stars than expected from a standard 
IMF. However, the discrepancy may also be attributed to different ages and star formation 
histories for the Cep 0B3b and ONC clusters. 

6. We show that the HD 217086 (07n) star with its known wind terminal velocity and 
mass-loss rate is unable to produce the bright soft X-ray diffuse emission seen in the field. 
An alternative explanation — the X-ray diffuse emission originates from the unresolved PMS 
point sources in the region — is suggested by (1) the spatial coincidence of the diffuse emission 
with the detected point source population, (2) the similarity of the plasma temperature of 
the diffuse emission with that of the combined weak PMS Orion-COUP sources, (3) the 
equivalence of the amount of the observed X-ray luminosity for the diffuse emission with 
that of the estimated undetected members of the unobscured population. 

7. We establish that the X-ray emission from the HD 217086 (07n) star is a constant 
soft-spectrum emission at the level log(Lx/-^6oO ~ —7 that is expected from the standard 
model involving many small shocks in an un-magnetized radiatively accelerated wind. 

8. From the evolution of an adaptively smoothed median energy during the unusually 
hot and strong flare in source # 294 we infer the slope of the trajectory in the temperature- 
density diagram, estimate the lower limit on the size of the flaring structure, and suggest 
that the source # 294 may join a number of similar super-flare sources seen in Orion clouds 
where long magnetic structures may connect the protoplanetory disk to the stellar surface. 
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Table 1. Basic Source Properties 



Source Position Extracted Counts Characteristics 



Seq 


CXOCEPB 


"J2000 


'5j2000 


e 


Net 


ANet 


Bkgd 


Net 


PSF 


Signif 


Anom* 


Var" 


EffExp 


Med E 


# 




(deg) 


(deg) 


(') 


Full 


Full 


Full 


Hard 


Frac 








(ks) 


(keV) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


1 


225531.25+624451.7 


343.880235 


62.747708 


10.22 


10.2 


4.0 


1.8 


2.9 


0.70 


2.2 


ge.. 




11.7 


1.70 


2 


226633.30+624457.2 


343.888775 


62.749237 


10.06 


9.1 


4.9 


8.9 


0.0 


0.90 


1.7 




a 


25.2 


1.16 


3 


225535.75+624235.4 


343.898960 


62.709839 


8.86 


8.0 


3.9 


3.0 


0.2 


0.90 


1.8 






21.2 


1.26 


4 


225537.53+624128.1 


343.906375 


62.691159 


8.38 


101.7 


10.8 


3.3 


67.9 


0.89 


9.0 




b 


26.3 


2.66 


5 


225538.33+624117.8 


343.909731 


62.688296 


8.26 


542.3 


23.9 


3.7 


178.4 


0.89 


22.2 




c 


26.6 


1.64 


6 


225540.99+624625.9 


343.920829 


62.773882 


10.18 


27.6 


6.7 


9.4 


15.2 


0.90 


3.8 




a 


25.3 


2.65 


7 


225541.56+623905.9 


343.923167 


62.651662 


7.83 


6.4 


3.5 


2.6 


0.0 


0.90 


1.5 




a 


27.0 


1.01 


8 


225541.74+624232.1 


343.923948 


62.708936 


8.18 


17.1 


5.0 


2.9 


2.3 


0.90 


3.1 




a 


24.5 


1.45 


9 


226642.16+623742.3 


343.925698 


62.628442 


8.03 


6.6 


3.7 


3.4 


0.0 


0.90 


1.6 




a 


24.9 


1.20 


10 


225544.19+624122.1 


343.934156 


62.689499 


7.61 


29.4 


6.2 


2.6 


8.4 


0.89 


4.4 




a 


26.9 


1.38 


11 


225544.21+624213.4 


343.934247 


62.703740 


7.82 


19.0 


5.1 


2.0 


2.7 


0.89 


3.3 




a 


26.8 


1.20 


12 


225545. ()8 + 6238,")4. 6 


343.937848 


62.648524 


7.45 


69.6 


9.0 


2.4 


20.4 


0.89 


7.3 




c 


27.2 


1.48 


13 


225547.50 + 624157.6 


343.947954 


62.699356 


7.38 


28.1 


6.0 


1.9 


0.8 


0.89 


4.3 




a 


27.0 


1.34 


14 


225547.59+623614.8 


343.948310 


62.604131 


7.99 


38.6 


7.0 


3.4 


9.2 


0.90 


5.1 


g... 




26.2 


1.23 


15 


225549.82+623421.9 


343.957606 


62.572759 


8.82 


15.2 


4.9 


3.8 


12.5 


0.90 


2.8 




a 


26.8 


2.99 



CO 



Note. — Column 1: X-ray source number. Column 2: lAU designation. Columns 3-4: Right ascension and declination for epoch J2000.0 in degrees. Column 5: OfF-axis 
angle. Columns 6-7: Estimated net counts from extracted area in total energy band (0.5 — 8.0) keV and their Ifj errors. Column 8: Estimated background counts from 
extracted area in total energy band. Column 9: Estimated net counts from extracted area in hard energy band (2.0 — 8.0) keV. Column 10: Fraction of the PSF at the fiducial 
energy of 1.497 keV enclosed within the extracted area. Note that a reduced PSF fraction (significantly below 90%) may indicate that the source is in a crowded region or 
on the edge of the field. Column 11: Source significance in sigma. Column 12. Source anomalies described by a set of flags, see the text note^ below for details. Column 
13: Variability characterization based on K-S statistics, see the text note^ below for details. Column 14: Source's effective exposure. Column 15: The background-corrected 
median photon energy in the total energy band. 

^Source anomalies: g = fractional time that source was on a detector (FRACEXPO from mkarf) is < 0.9 ; e = source on field edge; p = source piled up; s = source on 

readout streak. 

^Source variability: a — no evidence for variability; b — possibly variable; c = definitely variable. No test is performed for sources with fewer than 4 total full-band 
counts. No value is reported for sources in chip gaps or on field edges. 



Note. — 



The full table of 431 Chandra sources is available in the electronic edition of the Journal. 



Table 2. 2MASS Counterparts 



Seq 


CXOCEPB 


2MASS 


Off 


J 


H 


Ks 


Ph' 


Con'' 


^V,7K 




logibol.IR 


# 








mag 


mag 


mag 


Flag 


Flag 


mag 


Mq 


L0 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


1 


225531.25+624451.7 


22553069+6244482 


6.16 


14.76 


13.85 


13.53 


AAA 


000 


1.84 


0.15 


-0.86 


2 


225533.30+624457.2 


22553312+6244596 


2.72 


13.93 


13.00 


12.69 


AAA 


000 


1.39 


0.27 


-0.57 


3 


225535.75+624235.4 


22553545+6242343 


2.29 


12.54 


11.67 


11.38 


AEE 


000 


0.46 


0.68 


-0.10 


4 


225537.53+624128.1 


22553745+6241280 


0.57 


12.67 


11.25 


10.47 


AAA 


000 


6.15 


1.81 


0.74 


5 
6 


225538.33+624117.8 
225540.99+624625.9 


22553834+6241177 


0.14 


12.99 


11.96 


11.44 


AAA 


000 


1.62 


0.63 


-0.14 


7 


225541.56+623906.9 


22554154+6239030 


2.92 


12.06 


11.72 


11.66 


AAA 


cOO 








8 


225541.74+624232.1 


22554177+6242322 


0.19 


14.11 


13.17 


12.72 


AAA 


cOO 


1.49 


0.24 


-0.63 


9 


225542.16+623742.3 


22554205+6237426 


0.84 


15.28 


14.32 


14.01 


AAA 


000 


2.47 


0.12 


-1.01 


10 


225544.19+624122.1 


22554423+6241220 


0.29 


14.35 


13.45 


13.15 


AAA 


000 


1.51 


0.18 


-0.74 


11 


225544.21+624213.4 


22554407+6242130 


1.05 


13.94 


13.04 


12.80 


AAA 


000 


1.14 


0.25 


-0.61 


12 


225545.08+623864.6 


22554507+6238647 


0.07 


14.06 


13.13 


12.86 


AAA 


000 


1.45 


0.24 


-0.62 


13 


225547.50+624167.6 


22554749+6241676 


0.12 


12.96 


12.05 


11.78 


AAA 


000 


0.70 


0.53 


-0.25 


14 


225547.59+623614.8 


22554765+6236162 


0.68 


13.21 


12.29 


11.87 


EAA 


000 


0.69 


0.44 


-0.35 


15 


225549.82+623421.9 























00 



Note. — Column 1 : X-ray source number. Column 2: lAU designation. Column 3: 2MASS source. Column 4: C handr a-2M ASS 
positional offset. Columns 5-7: 2MASS JHKg magnitudes. Columns 8-9: 2MASS photometry quality, and confusion and contamination 
flags. See text notcs^''^ below for more details. Columns 10-12: Visual absorption, stellar mass, and bolometric luminosity derived from 
color-magnitude diagram by de-reddening 2MASS photometric colors to the 1 Myr PMS isochrone. 

^Photometry quality flag with three characters refer to the JHKs bands: A — very high significance detection (> 10 SNR,); B — high 
significance detection {> 7 SNR); C — moderate significance detection {> 5 SNR); D — low significance detection; E — point spread fitting 
poor; F — reliable photometric errors not available; U — upper limit on magnitude (source not detected); X — source is detected but no 
valid brightness estimate is available. 

'^2MASS confusion and contamination flag with three characters refer to the JHKs bands: = source unaffected by artifacts; b = 
bandmerge confusion (possible multiple source); c = photometric confusion from nearby star; d = diffraction spike confusion from nearby 
star; p = persistence contamination from nearby star; s — electronic stripe from nearby star. 



Note. — The full table of 431 Chandra sources is available in the electronic edition of the Journal. 



Table 3. X-ray Spectroscopy and Luminosities 



Seq 


CXOCBPB 


log Nh 


kT 




A log Lh 


logife.c 


Alogii,,c 


log Lt 


A log Lt 


log Lt,c 


A log Lt,c 


MMD 


# 






keV 


erg 


erg 


erg 


erg 


erg 


erg 


erg 


erg 


Flag 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


1 


225531. 25+624461. 7 


21.67 


4.70 


29.75 


0.34 


29.77 


0.33 


29.88 


0.18 


30.01 


0.26 


wlO 


2 


225533.30+624467.2 


21.44 


1.10 


28.50 


0.40 


28.53 


0.39 


29.09 


0.21 


29.36 


0.26 


wlO 


3 


225535.75+624235.4 


20.97 


1.90 


28.76 


0.51 


28.77 


0.51 


29.18 


0.21 


29.26 


0.20 


wlO 


4 


225537.53+624128.1 


22.14 


4.65 


30.52 


0.07 


30.58 


0.15 


30.67 


0.06 


30.83 


0.38 


blO 


B 


225538.33+624117.8 


21.69 


2.63 


30.90 


0.10 


30.93 


0.07 


31.11 


0.06 


31.29 


0.21 


blO 


6 


225540.99+624626.9 


22.61 


1.88 


29.71 


0.11 


29.92 


0.25 


29.76 


0.10 


30.38 


0.60 


baO 


7 


225541.56+623906.9 


19.00 


0.81 


27.89 


0.64 


27.89 


0.64 


29.06 


0.22 


29.06 


0.21 


vfO 


8 


225541.74+624232.1 


21.47 


2.10 


29.21 


0.33 


29.23 


0.32 


29.49 


0.17 


29.67 


0.22 


wlO 


9 


225542.16+623742.3 


21.69 


1.00 


28.37 


0.54 


28.43 


0.52 


28.94 


0.24 


29.39 


0.36 


wlO 


10 


225544.19+624122.1 


21.77 


0.69 


28.45 


0.21 


28.51 


0.19 


29.39 


0.13 


30.07 


0.31 


blO 


11 


225544.21 + 624213.4 


21.36 


1.40 


28.93 


0.25 


28.95 


0.24 


29.42 


0.15 


29.57 


0.17 


wlO 


12 


225545.08+623864.6 


21.30 


1.32 


29.36 


0.16 


29.38 


0.14 


29.89 


0.09 


30.06 


0.12 


blO 


13 


225547.50+624167.6 


21.24 


1.10 


28.69 


0.21 


28.70 


0.20 


29.38 


0.13 


29.64 


0.13 


blO 


14 


225547.59+623614.8 


21.02 


1.37 


29.26 


0.18 


29.27 


0.17 


29.82 


0.11 


29.91 


0.10 


blO 


16 


225549.82+623421.9 


22.21 


15.00 


29.74 


0.18 


29.80 


0.21 


29.77 


0.14 


29.96 


0.38 


vaO 



Note. — Column 1: X-ray source number. Column 2: lAU designation. Columns 3-4: Estimated column density and plasma energy. Columns 5-12: Observed 
and corrected for absorption X-ray luminosities, obtained from our spectral analysis (see section 4), and their errors estimated from simulations (see section 
4.3). h — hard band (2.t) — 8.0) kcV, t — total band (0.5 — 8.0) keV, c = corrected for absorption. Column 13: Three-part flag summarizing spectral method, 
source membership, and if-disk star information: 

b.. = spectral results from spectral fits of 174 relatively "bright" sources with > 20 net counts. 

w.. = spectral results from N^ jj^ — MedE interpolation for 212 "weak" sources with known N^ jj^. 

v.. = spectral results from fits of un-grouped spectra for 45 "weak" sources with unknown iVjy /jj. 

.e. = possible members of the embedded cluster. 

.1. — possible members of the lightly absorbed cluster. 

.o. — other possible members, found in the molecular cloud, not related to the embedded cluster, 
.a. = AGN candidates, 
.f. = field star candidates. 

..d = possible if-band excesses associated circumstellar disks. 
..0 = stars with no /f-band excesses or unknown. 



Note. — 



The full table of 431 Chandra sources is available in the electronic edition of the Journal. 
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Table 4. Power Law Fits for Extragalactic Candidates 



Seq 


CXOCEPB 


log Nh 


r 






# 




cm"*^ 


erj 


g cm 


-2 


(1) 


(2) 


(3) 


(4) 


(S) 




6 


226640.99+624625.9 


22.06 


1.71 


-13.73 




15 


225549.82+623421.9 


22.30 


1.06 


-13.84 




27 


225555.20+624356.6 


22.22 


2.02 


-13.01 




34 


225557.13+623751.7 


22.71 


-1.32 


-14.12 




41 


225601.45+623923.1 


22.11 


0.11 


-14.13 




45 


225603.10+623714.7 


23.11 


1.17 


-14.13 




64 


226608.82+623440.5 


20.00 


-0.11 


-14.27 




86 


226614.07+623552.7 


21.88 


1.56 


-14.26 




109 


225620.22+624321.8 


20.00 


0.68 


-14.36 




123 


225623. 1 1+623557. 1 


22.53 


2.93 


-14.08 




132 


225624.66+624548.7 


22.81 


2.16 


-13.70 




144 


225627.11+623442.8 


23.14 


9.50 


-14.32 




209 


225641.33+624259.4 


22.17 


1.48 


-14.10 




215 


226642.23+624016.3 


22.17 


0.21 


-14.20 




218 


226642.78+623858.8 


20.00 


0.82 


-14.61 




220 


225642.80+623211.9 


21.48 


1.41 


-14.30 




227 


225643.88+623832.4 


21.60 


1.46 


-14.55 




244 


225647.41+623057.9 


23.22 


0.95 


-13.41 




305 


225658.56+624434.8 


22.66 


3.01 


-14.39 




398 


226726.61+624628.1 


22.46 


3.11 


-14.37 




401 


226727.44+624516.5 


21.67 


1.07 


-14.29 




415 


226734.11+624251.1 


22.08 


0.88 


-13.93 




417 


226738.06+624238.2 


22.87 


4.20 


-14.21 




429 


226761.63+624813.0 


23.38 


9.46 


-13.86 




Note. 


— Column 1: X-ray source number. 


Column 2: 


lAU desi 


gna- 


tion. Columns 3-.5: Column density, photon 


index and 


observed 


flux 


in hard (2.0 — 8.0) keV energy band from spl 


ine-like fit 


of un-grouped 



spectrum. 
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Table 5. Integrated X-ray luminosities of detected PMS stars 



Population 


Ndet 














10^^ erg 


10^^ erg 


10^^ erg 


10^^ erg 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Cep OB3b unobscured 


321 


3.97 ± 0.10 


7.38 ± 0.89 


2.64 ± 0.11 


2.82 ± 0.14 


Cep B obscured 


64 


0.50 ± 0.03 


2.68 ± 1.49 


0.39 ± 0.03 


0.62 ± 0.06 



Note. — Column 1: Type of Cepheus PMS population. Column 2-6: Number of X-ray detected in 
this Chandra observation PMS stars and their total X-ray luminosities. 
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Table 6. EM linear regression 



log 


-.Lx =A- 


h S X lof 


5(M/Mq) 


Lx 


A 


AA 


B AS 


(1) 


(2) 


(3) 


(4) (5) 




29.95 


0.08 


2.26 0.13 


Lh.c 


29.97 


0.08 


2.23 0.13 




30.24 


0.07 


1.99 0.11 


Lt.c 


30.43 


0.07 


1.99 0.11 



Note. — Column 1: Type of Lx- Col- 
umn 2-3: Constant parameter of the lin- 
ear regression and its standard deviation. 
Columns 4-5: Slope of the linear regres- 
sion and its standard deviation. 



-40- 



Table 7. Integrated X-ray luminosities of undetected unobscured PMS stars 



Log Mass 


Bin 




T.Lt 


E it,c 






(NW-W) 


(NW-W) 


(NW-W) 








10^° erg 


10^° erg s-"^ 


(1) 




(2) 


(3) 


(4) 




-0.5 


32 - 50 


5.62 - 8.79 


8.71 - 13.61 




-0.7 


115 - 90 


8.08 - 6.33 


12.52 - 9.80 




-0.9 


165 - 151 


4.64 - 4.25 


7.18 - 6.58 




-1.1 


102 - 167 


1.15 - 1.88 


1.78 - 2.91 




E 


414 - 458 


19.49 - 21.24 


30.20 - 32.90 



Note. — Column 1: Mass bin in log. Column 2: Estimated 
number of X-ray undetected stars for both, subtractions of "NW" 
and "W" local backgrounds. Columns 3-4: Expected observed 
and corrected for absorption X-ray luminosities. 
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Table 8. X-ray properties of the diffuse component and PMS populations 



Object 


OR 


Nh 




kT2 


Lt 


it.c 




c/ks 


10^^ cm~^ 


keV 


keV 


10^^ erg 


10^^ erg 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Cepheus diffuse 


33.5 


0.49 ± 0.13 


0.78 ± 0.14 




2.21 ± 0.07 


6.92 ± 2.01 


Cepheus point sources 


570.0 


0.33 ± 0.23 


0.61 ± 0.23 


4.51 ± 0.62 


54.34 ± 4.41 


76.74 ± 19.91 


COUP 326 weak sources 


43.1 


0.34 ± 0.15 


0.77 ± 0.18 




1.93 ± 0.06 


4.84 ± 2.10 



Note. — Column 1: Defines object: either Cepheus X-ray diffuse emission, or 321 X-ray detected point source 
members of the Cep OB3b unobscurcd population, or 326 weak (log c < 29) members of the COUP unobscurcd 
population. Column 2: Net count rate. Column 3: Column density inferred from X S PEC fits, using wabs(apec) 
model for the diffuse emission and wabs(vapec-|-vapec) for the point source population. Columns 4-5: Plasma 
temperatures inferred from XSPEC fits. Columns 6-7: Observed and absorption corrected X-ray luminosities in 
(0.5 — 8.0) keV band, inferred from the X-ray spectral fits. 
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Fig. 1. — i?-band image covering 0.5° x 0.5° of the Cepheus B (Cep B) and Cep 0B3b 
neighborhood from the Digital Sky Survey. The Chandra 17' x 17' ACIS-I is outhned by the 
square. Cep B, the hottest component of the Cepheus Molecular Cloud, is at the bottom-left 
corner of the Chandra field. To the north and west lies Cep 0B3b, the youngest of two 
subgroups of the stellar Cep 0B3 association. The interface between Cep B and Cep 0B3 is 
delineated by the HII region, Sharpless 155 (S 155). The most massive and optically bright 
star in our field, HD 217086 (07n), is labelled. North is up, east is to the left. 




R.A.|J2pOO> 

Fig. 2. — The adaptively smoothed Cepheus B ACIS-I image in the full 0.5 — 8.0 keV band. 
Smoothing has been performed at the 2.5(7 level and gray scales are logarithmic. 
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Fig. 3. — Lightly absorbed Chandra stars which are not resolved in the 2MASS catalog. 
Each panel shows a 1' x 1' neighborhood (1 pixel = 1") around Chandra sources (marked by 
a red + and labelled with their source number in yellow) lying in close proximity to bright 
2MASS stars. Stars registered in 2MASS catalog are marked by cyan o. 
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H-KJiiiagJ 

Fig. 4. — NIR color-color diagram shows JHK colors of all (~ 2700) 2MASS sources (grey 
dots) inside the Chandra ACIS-I field, and 385 2MASS sources with possible Chandra coun- 
terparts (black dots). The solid and dot-dashed curves show sites of intrinsic JHKg colors 
of main sequence and giant stars respectively. The dashed lines are reddening vectors orig- 
inating at MO V (left line) and M6.5 V (right line), and marked at intervals of Ay — 2 
mag. 
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12 3 4 

J - H [mag] 



Fig. 5. — NIR color-magnitude diagram for the same set of sources as in Figure 4. Black dots 
show Chandra sources with 2MASS counterparts, and grey dots show all 2MASS sources in 
the field. Thick black lines give the unabsorbed ZAMS (dotted line) and 1 Myr isochrones 
(solid line), from which Ay ^ 10 mag reddening vectors are shown for various star masses. 
The remaining curves show various expected subpopulations of foreground and background 
stellar populations obtained from the Besangon simulations. Details are given in footnote 4. 
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Fig. 6. — X-ray source classification. The first five rows classify sources based on X-ray 
photometry and NIR counterpart-photometry information. The next three rows indicate 
classes of sources with different spectral analysis methods apphed. The last row shows 
classes of sources hkely unrelated to the Cepheus region derived in §§5.1 and 5.2. 
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Fig. 7. — Spectra of five lightly absorbed X-ray sources with 20 < NetFull < 23 counts 
exemplifying a common ambiguity in spectral fitting. Each row shows alternative models 
which successfully fit the data; the right-hand panel shows the fit we use with kT values 
typical for COUP PMS stars. 




Fig. 8. — Spectra of five heavily absorbed X-ray sources with NetFull > 20 counts exem- 
phfying stable spectral fitting over different grouping schemes. The finest grouping is to the 
left and the coarsest available grouping is to the right. As a final fit for highly absorbed 
sources we thus arbitrarily report a fit from one of the available groupings. 
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Fig. 9. — Plot of absorbing column densities obtained from X-ray spectral fits {NH_xray) 
against column densities obtained from NIR photometry {Nh_ir)- The sample consists of 
129 bright soft and 16 bright hard sources (see row 5 of Figure 6). Nhjr is derived from 
the visual absorption Ay using Nh = 2.0 X 10^^ Ay (Ryter 1996) where Ay is obtained 
by dereddening source colors to the 1 Myr PMS isochrone in the 2MASS color-magnitude 
diagram. The dashed hnes show \og{NH_Kray/ Nh_[r} of -0.3, 0, and 0.3. 
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Fig. 10. — Cumulative spectra for a sample of 12 consecutive sources of 177 weak soft sources 
with available Nh,ir information (row 7 in Figure 6). Upper panels show observed counts 
(step function) and model spectrum (smooth curve); lower panels show model residuals. 
NH,xray was Set equal to Nhjr- The source number and its net counts appear at the top 
left'. 
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Fig. 11. — Cumulative spectra for a sample of 12 consecutive sources of 35 weak hard sources 
with available Nh,ir information (row 7 in Figure 6). Upper panels show observed counts 
(step function) and model spectrum (smooth curve); lower panels show model residuals. 
NH,xray was Set equal to Nhjr- The source number and its net counts appear at the top 
left'. 
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Fig. 12. — Plots of formal errors in X-ray luminosities Alog(Lx), derived from our Monte 
Carlo simulations (see §4.3), as a function of source strength \og{N etFull / P S F) . The panels 
show results for four measures of source luminosity: observed hard band luminosities (L/^; 
2.0 — 8.0 keV), hard band luminosities corrected for absorption {Lh,c), observed total band 
luminosities (L^, 0.5 — 8.0 keV), and total band luminosities corrected for absorption (Lt^c)- 
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Fig. 13. — Spatial distribution of different classes of X-ray sources superposed on a grey-scale 
low-resolution image of the 17' x 17' Chandra field. In each panel, the embedded cluster is 
outlined by a black dashed circle of 3' radius. An outline of the HIT region and ionization 
front, seen in the optical DSS image is marked by the green contour, (a) Distribution of all 
83 hard {MedE > 2 keV) X-ray sources including NIR-identified (cyan x ) and un-identified 
(red X ) sources. The four ACIS-I CCD chips are labelled; the molecular cloud lies mostly on 
CCD 0. (b) Distribution of all 348 soft (MedE < 2 keV) X-ray sources including probable 
members of the region (cyan x) and possible foreground stars (red x). (c) iV//-stratified 
distribution of all probable cloud members; i.e., excluding the 24 (red x) probable AGNs 
in CCDs 1, 2, 3 from panel (a) and the 13 (red x) probable field stars from panel (b). 
Color code: yellow x have logA^^i^ < 21.3, red x have 21.3 < log Nh < 21.6, cyan x have 
21.6 < log A^H < 21.9, and blue x have logA^^ > 21.9. (d) Distribution of all 29 iT-band 
excess sources (blue x). HD 217086 07n star is marked by the red dot. 
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Fig. 14. — Evaluation of extragalactic source contamination using the X-ray source log N — 
logS diagram. A^(> S) gives the number of sources in the three CCDs outside of the 
molecular cloud. 5" is the source flux in the hard 2 — 8 keV band. The black sohd line shows 
the observed distribution for the 24 heavily absorbed sources without 2MASS counterparts. 
The dashed line shows the expected distribution for an unobscured extragalactic population 
(Moretti et al. 2003). The grey lines show examples of extragalactic source simulations, 
which include source spectral variations, large scale Galactic absorption, absorption local to 
Cepheus, and variations in the Chandra field X-ray background. 
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Fig. 15. — Expanded view of the foreground star locus area on the NIR color-magnitude 
diagram shown in Figure 5. Chandra sources are marked by black dots. The x symbols 
show simulated star locations of foreground stars based on the Besancon Galactic star model. 
Thirteen likely X-ray foreground stars are labelled and outlined by circles. 
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Fig. 16. — Positions of previously reported young members of the region, which fall within 
the ACIS-I field. Four VLA compact radio sources (□) and 3 red NIR objects (x) from 
(Testi et al. 1995), 27 ROSAT sources from (Naylor & Fabian 1999) (o), and 5 PMS with 
optical spectra from (Pozzo ct al. 2003) (+). Lines are the same as in Figure 13: the dashed 
circle outlines the embedded cluster of X-ray stars, and the solid contour outlines the HII 
region and ionization front seen in optical DSS image. 
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Fig. 17. — Comparison between XLFs of the ONC and Cepheus populations (§7.1). The 
upper sohd histograms show the COUP unobscured cool sample of 839 stars (Feigelson et al. 
2005); the middle dashed histogram shows the Cep 0B3b unobscured sample of 321 stars; 
the lower solid histogram shows the Cep B obscured sample of 64 stars. X-ray luminosities 
are the same as in Figure 12. Error bars indicate 68% confidence intervals (la) from Monte 
Carlo simulated distributions when X-ray luminosities are randomly drawn from Gaussian 
distributions with source's measured logL^ and source's estimated errors AlogLx- The 
broken red hues are added to aid the eye; they are based on the shape of the ONC XLFs 
and scaled downward to match the Cep 0B3b unobscured population. 
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Fig. 18. — Detection probability for a point source at a given observed total-band luminosity, 
as derived from our Monte-Carlo simulations of the Cep 0B3 ACIS observation. Three 
completeness curves correspond to three different detection thresholds: most representative 
to the current observation (solid), more conservative (dot-dashed), and less conservative 
(dashed). Three vertical lines indicate corresponding completeness limits: 29.5 (solid), 29.55 
(dot-dashed), and 29.4 (dashed). 
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Fig. 19. — Comparison between IMFs of the Orion Nebula Cluster and the Cepheus X-ray 
sources where stellar masses for both clusters are derived from NIR photometry (§7.2). The 

upper two histograms show the ONC derived from the X-ray COUP sample: the dashed is 
the 839-star COUP unobscured, and the solid shows the addition of 82 A'^-bright members 
un-detected in X-rays. The middle dashed histogram is for the 321-star X-ray unobscured 
Cep 0B3b population obtained here. The lower dotted histogram is for the 64-star obscured 
Cep B population. IMF models are shown as broken lines: dashed-dotted lines show three 
scaled versions of the ONC IMF (Muench et al. 2002), and grey hues show three scaled 
versions of Galactic field IMF (Kroupa 2002). 
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Fig. 20. — X-ray luminosity versus stellar mass. The grey dots denote the COUP unob- 
scured cool star sample of 839 stars. To be consistent with the current work, the plotted 
stellar masses for COUP sample have been determined from NIR color-magnitude diagram 
similar to stellar masses derived for Cepheus stars. Red and blue dots represent the Cepheus 
unabsorbed and embedded populations, respectively. X-ray luminosities are the same as in 
Figure 12. 
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Fig. 21. — The soft-band (0.5 — 2 keV) Chandra image after removal of point sources with 
holes smoothed over using an adaptive smoothing algorithm. The pentagon indicates the 
extraction region for the diffuse X-ray emission, and the white line outlines the area used 
for background subtraction in spectral analysis. Dashed lines show the location of the HII 
region with ionization front and the Cepheus obscured population, and the x marks the 
location of the HD 217086 07n star. 
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Fig. 22. — Background regions used in estimation of the Cep 0B3b unobscured cluster KLF. 
This is a 3° X 3° image from the i?-band Digital Sky Survey. The ACIS-I field is outlined by 
the white square. The local number of 2MASS sources in a 127 arcmin^ region was measured 
at each of the circled locations. The two areas which we use for the background subtraction, 
with 814 sources in the Western field and 830 sources in the Northwestern field, are marked 
by white crosses. 







1 2 
J-H [mag] 



3 



Fig. 23. — Comparison of 2MASS J — H color distributions in the cluster, molecular cloud, 
and background fields. The solid hne shows PMS stars of the Cep 0B3b unobscured pop- 
ulation in the ACIS-I field. The dotted line shows absorbed 2MASS stars in the molecular 
cloud (region around 22'*58"*, -|-62°30' in Figure 22). The upper (lower) dashed lines show 
2MASS sources in the Western (Northwestern) background fields. 
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Fig. 24. — Comparison of IMFs assuming an age of 1 Myr (a) and 3 Myr (b). The dashed 
histogram shows the IMF of the Cep 0B3b unobscured X-ray population with masses derived 
from the NIR color-magnitude diagram (same as in Figure 19). The solid lines represent the 
IMF trends derived from the 2MASS population within the Chandra field after subtraction 
of KLFs from the Western (A) and Northwestern (o) background fields and conversion them 
to IMFs (§7.3). The dashed-dotted broken line and the grey line show scaled versions of the 
ONC IMF (Muench et al. 2002) and Galactic field IMF (Kroupa 2002) respectively 
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Fig. 25. — Cepheus composite X-ray spectra, (a) Spectrum of the diffuse component detected 
in the 127 arcmin^ region shown in Figure 21 (points with error bars) with spectral model 
(histogram), (b) The composite spectrum and spectral model for 321 PMS stars of the 
unobscured X-ray population. The lower plots show residuals between the data and models 
in units of x^- 




Fig. 26. — Composite X-ray spectra of the cool unobscured ONC population, (a) Subsample 
of 326 weak (log Lj^c < 29) sources, (b) Full sample of 820 sources with no pileup. 
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Fig. 27.— X-ray properties of the 07n star HD 217086 {Chandra source #240): (a) ACIS 
spectrum, and (b) lightcurve. 
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Fig. 28. — X-ray properties of the superfiare in Chandra source #294. (a) Lightcurve 
histograms of the source #294 with a bin size of a half of an hour for a full energy band 
0.5 — 8.0 keV (black histogram), soft band 0.5 — 2.0 keV (red histogram), and hard band 
2.0 — 8.0 keV (blue histogram). Black and blue curves represent unbinned but adaptively 
smoothed versions of the lightcurve with a window size of at least 256 and 570 counts 
respectively. Vertical labelled yellow lines give two characteristic decay ranges [s,c] (15.5 — 
19.9 ks) and [c,e] (19.9 — 27.8 ks), found from our further analysis of the adaptively smoothed 
median energies (Figure 31). (b) Spectra of four phases of the lightcurve with 1-T optically 
thin thermal plasma fits, subject to logA^^^ ~ 21.8 cm~^ absorption: quiescent level (cyan, 
kT ~ 1.7 keV), flare rise and peak (magenta, kT ~ 64 keV), early decay phase (green, 
kT ~ 12 keV), and late decay phase (black, kT ~ 4.5 keV). 
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Fig. 29. — The lightcurve histogram (black dots with uncertainties) with a bin size of 1.06 
hours for the COUP source # 1343, used here to vahdate our method of modeUing the 
Cepheus source ^^294 flare. Grey and black lines present adaptively smoothed versions of 
the lightcurve with a window size of at least 100 and 625 counts respectively. Vertical labelled 
dashed lines give the time range [s,e] (230 — 330 ks) over which our analysis of the adaptively 
smoothed median energies was performed. 
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Fig. 30. — X-ray spectral modelling of the superflare in ONC COUP source #1343, using 
the method of the adaptively smoothed median energy with a window size of at least 625 
counts over the time range [s,e] (230 — 330 ks), including the peak of the flare and the decay 
phase. The top left panel shows the evolution of the observed median energy (88 points), 
the top right panel - the evolution of the plasma temperature, the bottom left panel - the 
evolution of the emission measure. The bottom right panel shows the evolution of the flare 
in the log T — log y^E'M) plane with the best fitting decay as a solid line of the slope 
C = 1.86 ±0.29. 
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Fig. 31. — X-ray spectral modelling of the superflare in source #294, using the adaptively 
smoothed median energy with a window size of at least 570 counts over the time range [s,e] 
(15.5 — 27.8 ks), including the decay phase of the flare (see Figure 28). The top left panel 
shows the evolution of the smoothed median energy for the window sizes of 256 (grey) and 
570 (black) counts respectively. The top right panel is the evolution of the temperature (403 
points), the bottom left panel is the evolution of the emission measure. The bottom right 
panel shows the evolution of the flare in the log T — log \/{EM) plane with the two best 
fitting decays: the slope of C = 8.08 ± 0.96 over the [s,c] (15.5 — 19.9 ks) time range; and the 
slope of C = 1.40 ± 0.13 over the [c,e] (19.9 - 27.8 ks) range. 



